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The static vacuum Einstein equations are the equations 



o 

o , 

O . ur = D'^u, (0.1) 



Au = 0, 

on a Riemannian 3-manifold {M,g), with u a positive function on M. Here r denotes the Ricci 
curvature, the Hessian, and A = trD^ the Laplacian on {M,g). Solutions of these equations 
5^ ! define a Ricci- flat 4-manifold N, of the form N = M Xu oi N = M XuM., with Riemannian or 

Lorentzian metric of the form 



gN = 9M±u^dt\ (0.2) 

■ These equations are the simplest equations for Ricci-flat 4-manifolds. They have been extensively 

studied in the physics literature on classical relativity, where the solutions represent space-times 
outside regions of matter which are translation and reflection invariant in the time direction t. 
, However, with the exception of some notable instances, (c.f. Theorem 0.1 below), many of the 

global properties of solutions have not been rigorously examined, either from mathematical or 
^ physical points of view, c.f. [Br] for example. 

qh| This paper is also motivated by the fact that solutions of the static vacuum equations arise in the 

study of degenerations of Yamabe metrics (or metrics of constant scalar curvature) on 3-manifolds, 
bX). c.f. [Al]. Because of this and other related applications of these equations to the geometry of 

3-manifolds, we are interested in general mathematical aspects of the equations and their solutions 
which might not be physically relevant; for example, we allow solutions with negative mass. 

In this paper, we will be mostly concerned with the geometry of the 3-manifold solutions (M, g, u) 
' of (0.1), (i.e. the space-like hypersurfaces), and not with the 4-manifold metric. Thus, the choice of 

Riemannian or Lorentzian geometry on in (0.2) will play no role. This considerably simplifies the 
discussion of singularities and boundary structure, but still allows for a large variety of behaviors; 
c.f. [ES] for a survey on singularities of space-times. 

Obviously, there are no non-flat solutions to (0.1) on closed manifolds, and so it will be assumed 
that M is an open, connected oriented 3-manifold. Let M be the metric (or Cauchy) completion 
of M and dM the metric boundary, so that M = M U dM is complete as a metric space. 

In order to avoid trivial ambiguities, we will only consider maximal solutions of the equations 
(0.1). For example any domain Q in with the flat metric, and u a positive constant, satisfies 
(0.1). In this case, the metric boundary is artificial, and has no intrinsic relation with the 
geometry of the solution. The solution obviously extends to a larger domain, i.e. M^. Thus, we 
only consider maximal solutions {M,g,u), in the sense that {M,g,u) does not extend to a larger 
domain (M',g') D {M,g) with u > on M' . It follows that at the metric boundary dM of M, 
either the metric or u degenerates in some way or u approaches in some way, (or a combination 
of such). 



Partially supported by NSF Grants DMS 9505744 and 9802722. 

1 



2 



A classical result of Liclinerowicz [LI, p. 137] implies that if the metric {M,g) satisfying (0.1) is 
complete, (i.e. dM = 0), and « ^ 1 at infinity, then u = I and M is flat, i.e. or a quotient of R^. 
More generally, it is proved in [Al, Thm.3.2] that if {M,g) is a complete solution to (0.1), (hence 
u > everywhere), then {M,g) is fiat and u is constant, i.e. the assumption on the asymptotic 
behavior of u is not necessary, c.f. also Theorem 1.1 below. 

Thus, there are no complete non-trivial solutions to (0.1) and hence dM must be non-empty. 
The set formally given by 

S = {u = 0} C M, 

is called the horizon. It is closely related to the notion of event horizon in general relativity. More 
precisely, T, may be defined as the set of limit points of Cauchy sequences on (M, g) on which u 
converges to 0. Although dM 7^ 0, it is possible that S = 0. However, most solutions of physical 
interest do have S 7^ 0. In the framework of classical relativity, the non-triviality of a static vacuum 
solution, i.e. the non- vanishing of its curvature, is due to the presence of matter or field sources at 
dM, or 'inside' the horizon E in case {M,g) extends as a vacuum solution past E. 

It is natural to consider the situation where {M,g) is not complete and for which the metric 
boundary dM of M coincides with the horizon S. More precisely, we will say that {M,g) is 
complete away from the horizon S if for any sequence pi ^ p dM in the metric topology on 
M = M U dM one has u{pi) — > 0. Conversely if {pi} is a bounded sequence in M with u{pi) — > 
0, then the definition of {M,g,u) implies that a subsequence of {pi} converges to a point p G dM. 
Thus, dM = S is given by the Hausdorff limit of the s -levels -L^ of it as e ^ 0. 

While most solutions {M,g) of physical interest are complete away from S, there are many 
solutions for which this is not the case, c.f. §2 for further discussion. In such examples, the 
curvature typically blows up within a finite distance to S. Among the solutions which are complete 
away from E, most all are singular at E, again in the sense that the curvature of the metric g 
blows up on approach to S. This is closely related to the fact that the equations (0.1) are formally 
degenerate at S. Thus, in general, even when M is complete away from E, the metric completion 
M U S need not be a smooth manifold with boundary. 

We will say that (M,g,u) extends smoothly to S, if (i): the set S is a smooth surface and the 
partial completion M' = MUl]cMisa smooth manifold with boundary, (ii): the metric g extends 
smoothly to a smooth Riemannian metric on M', and (iii): the potential u extends smoothly to S 
with E = {n = 0}; smoothness here means at least C^. Note that it might be possible that E has 
infinitely many components, or non-compact components of infinite topological type. In any case, 
one immediate strong consequence of (0.1) is that if g extends smoothly to S, then at S one has 

D^u\j: = 0. (0.3) 

Thus, S is a totally geodesic surface in M and |Vit| is a non-zero constant on each component of S, 
c.f. Remark 1.5 below. Observe that if {M,g) is smooth up to E, and complete away from S, then 
the isometric double of M across S is a smooth complete Riemannian manifold. The harmonic 
function u extends smoothly across S as harmonic function, odd w.r.t. reflection in E. 

By far the most significant solution of the static vacuum equations is the Schwarzschild metric, 
of mass m, given on the space-like hypersurface M by 

gs = {l- ^)-^dr^ + r^dsl„u = (1 - ^)V2,^ > 2m. (0.4) 

This metric models the vacuum exterior region of an isolated static star or black hole. It is 
a spherically symmetric metric on M = (2m, 00) x 5*^ and has E given by a (totally geodesic) 
symmetric S^, of radius 2m. The mass m is usually assumed to be positive, but we will not make 
this assumption here. Thus, we allow m < 0. Of course if m = 0, then gs is just the fiat metric 
with u= 1. If m < 0, (0.4) is understood to be defined for r > 0. 
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The Schwarzschild metric is asymptotically flat in the sense that there is a compact set K in 
M such that M \ K is difFeomorphic to \ B{R), and in a suitable chart on M \K, the metric 
approaches the Euclidean metric at a rate of 1/r, i.e. 

gij = il + ^)Sij + 0{l/r% (0.5) 

with curvature decay of order l/r^,r = \x\, and with m G M. The function u (up to a multiplicative 
constant) has the asymptotic form u = 1 — ^ + 0(r~^) with \Vu\ = 0(l/r^). A triple {M,g,u) 
satisfying these conditions is called asymptotically flat. 

We note the following remarkable characterization of the Schwarzschild metric. 

Theorem 0.1. (Black-hole uniqueness), [ll],[Ro],[BM] 

Let (M, g, u) be a solution of the static vacuum Einstein equations, which is smooth up to E 
and complete away from S. // S is a compact, (possibly disconnected) surface and {M,g,u) is 
asymptotically flat, then (M, g, u) is the Schwarzschild metric, with m > 0. 

The hypothesis that the space-like hypersurface (M, g, u) is asymptotically flat is very common 
in physics. Namely, in modeling a static space-time outside an isolated, i.e. compact, field or 
matter source, it is natural to assume that in regions far away from the source the geometry of 
space approximates that of R^, i.e. empty space. 

Nevertheless, mathematically the asymptotically flat assumption is quite strong in that it severely 
restricts both the topology and geometry of (M, g) outside a large compact set. Further, the physical 
reasoning above presupposes that there are no complete non-flat solutions to the vacuum equations 
(0.1), i.e. that a static gravitational field is non-empty solely due the presence of matter somewhere. 

Remark 0.2. This latter issue in fact led Einstein to hypothesize that space AI is compact, in order 
to avoid dealing with 'artificial' boundary value problems at infinity, c.f. [E, p.98ff.]. This issue, 
closely related to Mach's Principle, is discussed in some detail in work of Lichnerowicz, (Propositions 
A and B of [LI, §31] and [L2, Ch.II]); c.f. also [MTW, §21.12,],[Ri, §9.12] for instance. 

As remarked above, there are in fact no complete non-trivial static vacuum solutions, so that the 
asymptotically flat assumption in Theorem 0.1 may not be unreasonable, c.f. however [El], [G]. In 
fact, the main result of this paper is that this hypothesis is not necessary in most circumstances; 

it follows from much weaker assumptions. 

To explain this, we first need to consider a weakening of the condition that dM is compact. Let 

t{x) = dist{x,dM). (0.6) 

Define dM to be pseudo- compact if there is a tubular neighborhood U of dM whose boundary 
dU in M, dU fl M, is compact, i.e. {t{x) = So} is compact, for some Sg > 0, (and hence all 
< So < oo). As will be seen in §2, there are numerous examples of static vacuum solutions with 
dM pseudo-compact but not compact. 

Let E be an end, i.e. an unbounded component, oi M \ U. The mass of E may be defined by 

niE = lims-^oofnE^s) = lims-*oo-r~ I < Vlogu, Vt > dA, (0-7) 

Jse(s) 

where Se{s) = t~^{s) fl E. Since u is harmonic, log u is superharmonic, so that the divergence 
theorem implies that mE{s) is monotone decreasing in s. Hence the limit (0.7) is well-defined, 
(possibly — oo). Note that the static vacuum equations are invariant under multiplication of u by 
positive constants. We use the logu term in (0.7) in place of u so that the mass is independent of 
this rescaling in u. 

The following is the main result of this paper. 
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Theorem 0.3. Let {M, g, u) be a solution of the static vacuum Einstein equations. 

(i) . Suppose dM is pseudo- compact. 

Then M\U has only finitely many ends {Ei\, < i < q < oo. Supposing q > and (i) holds, 
let E G {Ei} he any end of M satisfying 

(ii) . u{xi) > Uo, for some constant > and some sequence Xi (z E with t{xi) oo. 

Then the end {E, g, u) is either asymptotically flat, or small, in the sense that the area growth 
of geodesic spheres satisfies 



Further, if niE 7^ and supEU < 00, then the end E is asymptotically flat. 

We make several remarks on this result. First, if dM is not pseudo-compact, then it is easy to 
construct static vacuum solutions for which M\U has infinitely many ends, c.f. the end of §1 
or §2(1). Further, there are examples of ends E with compact boundary, on which (ii) does not 
hold and which are neither asymptotically flat nor small in the sense of (0.8), c.f. Example 2.11. 
Thus, both hypotheses (i) and (ii) are necessary in Theorem 0.3. There are also examples with M 
compact, c.f. §2(1), so that M\U may have no ends, {q = Q). 

On the other hand, both alternatives in Theorem 0.3, namely asymptotically flat or small ends, 
do occur. Asymptotically flat ends satisfy areaSE{s) ~ 47rs^ as s — > 00, while small ends have 
small area growth. For example, (0.8) implies, at least, that there is a sequence Sj — >■ 00 such that 



for any flxed e > 0. All known examples of solutions satisfying (0.8) are topologically of the form 
(M^ \ B)x outside a compact set and the geodesic spheres have at most linear area growth. The 
main example of a static vacuum solution with a small end is the family of Kasner metrics, c.f. 
Example 2.11 below. 

To illustrate the sharpness of the last statement in Theorem 0.3, we construct in Remark 3.8 a 
(dipole-type) static vacuum solution (M, g^ u) with a single end E on which tue = 0, supE < 00, 
and which satisfies (0.8). Similarly, Example 2.11 provides static vacuum solutions with tue > 0, 
supEU = 00 and satisfying (0.8). Thus, the last result is also sharp. 

Assuming dM is pseudo-compact, it is easy to see that if > then (ii) holds, so that (ii) may 
be replaced by the assumption > 0. Thus, for the physically very reasonable class of solutions 
such that dM is pseudo-compact, niE > for all ends E, and u is bounded, all ends E of M are 
asymptotically flat. We also point out that if M is complete away from S, then (ii) holds at least 
on some end E. 

The proof of Theorem 0.3 gives some further information on the asymptotic structure of the 
small ends. For instance, the curvature in such ends decays at least quadratically, c.f. (1.3), and 
in the geodesic annuli AE(isi,2si) in E, the metric approaches in a natural sense that of a Weyl 
solution, i.e. a static axisymmetric solution, as Si — 00. Thus, asymptotically, small ends have at 
least one non-trivial Killing field. However, it is not known for example if the metric is asymptotic 
to a unique Weyl solution, or even if small ends are necessarily of finite topological type. 

It would also be of interest to prove that (M, g, u) has a unique end if dM is pseudo-compact. 
However, we have not been able to do this without further assumptions, c.f. Remark 3.10. 

The definition of asymptotically fiat rules out the possibility that M is smooth up to S and 
complete away from S, with S non-compact, since for example u = on S while u — > 1 on any 
asymptotically flat end. Consider however the following (Bl) solution. 




(0.8) 



areaSE{si) < Sj • {logsi) 



gsi = (1 - ^)-^dr2 + (1 - ^)d(f>'' + r^d^^ 



(0.9) 
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where 0/4m G [0, 27r],^ G [0, 7r],r > 2m > 0, with u = r ■ sin{9). Note that the potential u is 
unbounded. 

The Bl metric is just a 3-dimensional sHce of the 4-dimensional Schwarzschild metric {N,gs + 
u^dO"^), {t changed to 9 in (0.2)), obtained by dividing N by an C Isoni{S'^) orthogonal to dO; 
this is the slice 'orthogonal' to the usual slice giving the Schwarzschild metric (0.4). This metric has 
S given by two disjoint, isometric copies of B?, each of positive Gauss curvature and asymptotic to 
a fiat cylinder. It is smooth up to S and complete away from S. The metric is globally asymptotic 
to the flat metric on x S^, again with curvature decay of order 0(l/r^) and with u of linear 
growth in distance to S. Such solutions will be called asymptotically cylindrical. In fact a large 
class of Weyl solutions have 'dual' solutions in this sense which are asymptotically cylindrical, c.f. 
Remark 2.9. 

This paper is organized as follows. Following discussion of some general topics on static space- 
times in §1, we analyse in some detail the class of Weyl vacuum solutions in §2. Several new 
results on the structure of these solutions are given; for instance Proposition 2.2 gives a new 
characterization of Weyl metrics. In addition, some efforts have been made to give a reasonably 
clear and organized account of the breadth of possibilities and behavior of Weyl metrics, since their 
treatment in the literature is rather sketchy and since they serve as a large class of models on which 
to test Theorem 0.3. Theorem 0.3 is proved in §3, and the paper concludes with several remarks 
on generalizations, and some open questions. 

I would like to thank the referee for suggesting a number of improvements in the exposition of 
the paper. 

1. Background Discussion. 

Let {M,g,u) be an open, connected oriented Riemannian 3-manifold and N = Af x„ R or 
= M x „ S^, as in (0.2). Thus N represents a static space-time and M C A'" is totally geodesic. 
The Einstein field equations on are 

rN-^-gN = T, (1.1) 

where T is the (M or S'-'^-invariant) stress-energy tensor. (We are ignoring physical constants here). 
These equations may be expressed on the space-like hypersurface M as the system 

r -u-^D'^u + {u~^Au- Is) ■ g = Th, (1.2) 
i« = Tv, 

where Th is the horizontal or space-like part of T and Ty is the vertical or time-like part of T. 
These are the equations on M for a Lorentzian space-time A'^; in case A^ is Riemannian, the first 
equation is the same while the second is -^s = Ty- When T = 0, one obtains the vacuum equations 
(0.1), which are the same for Lorentzian or Riemannian signature. A common example, with T ^ 
is a static perfect fluid, c.f. [Wd, Ch.4], with T given by 

T = (/X + p)dt'^ + pg, 

where /i, p are time independent scalar fields representing the energy density and pressure respec- 
tively. The equations (1.2) imply that the full Riemann curvature i?jv of A^ is determined by r, u 
and T. 

The horizon S = {-u = 0} corresponds formally to the fixed point set of the action on N 
and requires special consideration. For example, even if M is smooth up to E the Riemannian 4- 
manifold {N,gN) might not be smooth across S, even though the curvature Rn is smooth. Namely, 
assuming the parameter t G [0, 27r), if | Vu|s ^ 1, then N has cone singularities (with constant 
angle by (0.3)) along and normal to the totally geodesic submanifold H C N. (This issue does not 
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arise for Lorentzian metrics). By multiplying the potential fTinction it by a suitable constant, one 
can make the metric smooth across any given component of S; one cannot expect however in 
general that this can be done simultaneously for all components of E, if there are more than one. 
This issue will reappear in §2. 

The following result is proved in [Al, Cor. A. 3]. It implies, (by letting t co), that if {M,g,u) 
is a complete solution to the static vacuum equations with u > everywhere, then M is flat, and 
u is constant. 

Theorem 1.1. Let (M, g, u) he a solution to static vacuum equations (0.1). Let t(x) = dist(x, 
dM) as in (0.6). Then there is a constant K < oo, independent of (M, g, u), such that 

H(x)<j^,|Vi»9u|(x)<^. (1.3) 

Remark 1.2. The same result has recently been proved for stationary vacuum solutions, i.e. space- 
times admitting a complete time-like Killing field, c.f. [A2] . We will discuss elsewhere to what extent 
Theorem 0.3 generalizes to stationary vacuum solutions. 

Recall from §0 that the potential function u of a static vacuum solution may be freely renormal- 
ized by arbitrary positive constants; hence the appearance of log u in (1.3), as in (0.7). 

Theorem 1.1 implies that the curvature of {M,g) is controlled away from dM, and hence the 
local geometry of solutions is controlled away from dM by lower bounds on the local volume or 
injectivity radius. More precisely, we have the following results which are essentially a standard 
application of the Cheeger-Gromov theory of convergence/collapse of Riemannian manifolds, c.f. 
[CG], [A3, §2], or also [P, Ch.l2] for an introduction to these results. Further details of the proofs 
of these results are given in [A2] , (for the more general class of stationary space-times) , and also in 
[Al,App.]. 

Lemma 1.3. (Non-Collapse). Let {Mi,gi,Ui) be a sequence of solutions to the static vacuum 
equations (0.1). Suppose 

\ri\ < A,diamMi < D,volMi > fo, 

and 

dist{xi, dMi) > 5, 

for some Xi G Mj and positive constants Vo, D, S. Assume also that Ui is normalized so that 

Ui{Xi) = 1. 

Then, for any e > sufficiently small, there are domains Ui C Mj, with e/2 < dist{dUi, dMi) < e, 
and Xi G Ui such that a subsequence of the Riemannian manifolds {Ui,gi) converges, in the C°° 
topology, modulo diJJeom,orphisms, to a limit manifold (U, g), with limit function u and base point 
X = lim Xi. The triple (U, g, u) is a smooth solution of the static equations with u(x) = 1. 

Lemma 1.4. (Collapse). Let {Mi,gi,Ui) be a sequence of solutions to the static vacuum equations 
(0.1). Suppose 

\ri\ < A, diamMi < D, volMi 

and 

dist{xi, dMi) > 6, 

for some Xi G Mj and constants A, D, S. Assume also that Ui is normalized so that Ui{xi) = 1. 

Then, for any e > sufficiently small, there are domains Ui C Mj, with e/2 < dist{dUi, dMi) < e 
with Xi € Ui, such that U is either a Seifert fibered space or a torus bundle over an interval. In 
both cases, the gi- diameter of any fiber F, (necessarily a circle or torus T^), goes to as i ^ oo, 
and 7ri(F) injects in TTiiUi). 
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Consequently, there are infinite Z or Z0Z covers Ui ofUi, such that {Ui,gi,Xi} does not collapse 
and hence has a subsequence converging smoothly to a limit (U, g, x) of the static vacuum equations 
with X = Urn x\, x\ a lift of Xi to Ui. The limit {U, g, x) admits a free isometric M or M © M action, 
(c.f. %2), which also leaves the potential function u invariant, and u{x) = 1. 

In studying static solutions, it is often very useful to consider the conformally equivalent metric 
g = ■ g on M. An easy calculation using the behavior of Ricci curvature under conformal 
deformations, c.f. [Bes, p. 59], shows that the Ricci curvature f of g, in the vacuum case (0.1), is 
given by 

f = 2{dloguf > 0. (1.4) 
Further, if A denotes the Laplacian of g, then 

Klogu = 0. (1.5) 

The equations (1.4)-(1.5) are equivalent to the static vacuum equations (0.1). Since these equations 
are invariant under the substitution u — ^ — u, it follows that if (M, g, u) is a static vacuum solution, 
then so is {M,g',u~^), with g' given by 

/ 4 
9 = u -g. 

Similarly, observe that if (iV, ^jv) is the associated Ricci-flat 4-manifold (0.2), then 

ANlogu = 0. (1.6) 

Here and below, log always denotes the natural logarithm. 

We discuss briefly some of the simplest static vacuum solutions: 

Levi-Civita Solutions. 

There are 7 classes of so-called degenerate static vacuum solutions, where the eigenvalues Aj of 
the Ricci curvature r satisfy Ai = A2 = -2A3, called A1-A3, B1-B3, C, c.f. [EK, §2-3.6]. The 
B metrics are dual to the A metrics, as mentioned in §0, c.f. §2 for details. The Al metric is 
the Schwarzschild metric. It is of interest to examine the A2 metric, given in standard cylindrical 
coordinates on by 

gA2 = z'^idr^ + {sinh^r)d(t>^) + - l)~^(^2^ (1-7) 

with u = (^ — 1)^/^ and z G [0, 2m], m > 0. The horizon S = {u = 0} is given by the set {z = 2m} 
and hence is the complete hyperbolic metric H^, with curvature — (2m)~^. It is easily verified that 
the A2 metric is smooth up to S. However, S 7^ dM; the set {z = 0} is at finite distance to S, and 
so dM has another (singular) component obtained by crushing (compact subsets of) the hyperbolic 
metric to a point. 

Let r be any discrete group of hyperbolic isometries. Then T extends in an obvious way to a 

group of isometries of (7^12 • The uniformization theorem for surfaces implies that any orientable 
surface except 5^ and T^, including surfaces of infinite topological type and infinitely many ends, 
admits a complete hyperbolic metric, i.e. is the quotient H^/T, for some F. Hence, any such surface 
and hyperbolic metric can be realized as the horizon S of a static vacuum solution. 

Topologically, for S = H'^/T, we have M = S x /. Hence for example if S has infinitely many 
ends, then M also has infinitely many ends; in particular, this shows that the hypothesis that dM 
is pseudo-compact in Theorem 0.3 is necessary. 

The A3 metric is 



£,^3 = z^{dr^ + r^dcf)^) + zdz^ , 



(1.8) 
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with u = z~^l'^ > 0, r > 0. Hence E is empty in this case - it occurs at infinity in the metric. This 
metric may be reahzed as a pointed hmit of the A2 metric as m — > cx) and also as a hmit, in a 
certain sense, of the Al metric, c.f. Example 2.11. 

Remark 1.5. The discussion above raises the natural question if any orientable connected Rie- 
mannian surface can be realized as the horizon of a static vacuum solution, smooth up to S, 

which is defined at least in a neighborhood of S. In general, this appears to be unknown. 

Observe that any complete constant curvature metric on an orientable surface can be realized in 
this way. The Schwarzschild metric gives the constant curvature metric on 5^, the quotients of the 
A2 metric give all hyperbolic surfaces, and quotients of the flat metric, with u a linear function, 
give all flat metrics on a surface, (T^,^^ x M.or M^). 

Geroch-Hartle show in [GH] that any rotation-invariant metric on S'^ or can be realized at 
the horizon. Except for the Schwarzschild metric, such solutions are not complete away from S. 

Observe that the full 1-jet of (M, g) at S (assumed connected) is determined solely by the surface 
metric (S, (7), since S is totally geodesic and, renormalizing u if necessary, |Vii| = 1 on S by (0.3). 
Observe also that one cannot have |Vit| = on E, since u is harmonic and the divergence theorem 
applied to a small neighborhood of S would imply that u = on [/, which is ruled out. 

On the other hand, the metric {M,g) is not uniquely determined by its boundary values {T,,g). 
Namely, the flat metric on is realized by the flat vacuum solution M = T^xR+,u = t = dist(T,, •) 
and also (locally) by a non-flat metric, c.f. [T],[P]. Similar remarks hold for local perturbations of 
the Schwarzschild metric, c.f. [GH]. 

2. Weyl Solutions. 

A large and very interesting class of explicit solutions of the static vacuum equations are given by 
the Weyl solutions [W], c.f. also [EK,§2.3-9] or [Kr, Ch. 16-18]. In fact, it appears that essentially 
all known explicit solutions of the static vacuum equations are of this form. Since the literature 
on these solutions is not very organized or rigorous, especially regarding their global structure, 
we discuss these solutions in some detail. These metrics will also illustrate the necessity of the 
hypotheses in Theorem 0.3. 

Definition 2.1. A Weyl solution is a solution {M,g, u) of the static vacuum equations (0.1) which 
admits an isometric 'R-action, i.e. a non-zero homomorphism M — > Isom{M), leaving u invariant. 

Apriori, the topology of a Weyl solution could be quite non-trivial; for example M could be 
any Seifert fibered space. The first result shows that only the simplest topology (and geometry) 
is possible. For the moment, we exclude any possible fixed point set of the R-action from the 
discussion. 

Proposition 2.2. Let (M, g, u) be a Weyl solution with R-action without fixed points, which does 
not admit a (local) free isometric M x M action. Then the universal cover {M,g) of {M,g) is a 
warped product of the form 

M = VXfR,g = gv + f''d4>'', (2.1) 

with (V, gv) a Riemannian surface and f a positive function on V. The R-action on M is by 
translation on the second factor. 

Proof: This result, whose proof is purely local, is a strengthening in this situation of a well- 
known result in general relativity, Papapetrou's theorem, c.f. [Wd,Thm.7.1.1], which requires 
certain global assumptions, (e.g. smoothness up to S). 

Let K denote the (complete) Killing field generated by the R-action on M, and / = \K\. We 
may assume that u is not a constant function on M, since if u is constant, the metric is fiat, and so 
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admits a local R x M action. Since (M, g, u) is real-analytic, u is not constant on any open set in M. 
We thus choose a neighborhood U of any point p where Vu{p) 7^ on which \ Vu\ > 0. Define ei by 
ei = Vit/|Vit|, and extend it to a local orthonormal frame 61,62,63 for which 63 = K/\K\ = K/f. 
Note that this is possible since u is required to be invariant under the flow of if, so that 

<Vu,K>=0. (2.2) 

In U, the metric g may be written as 

g = 7r*gv + f\dcl> + ef, (2.3) 

where tt : U V is a Riemannian submersion onto a local surface {V, gv), 6 is a connection 1-form, 
K = d/d(j) and / is a function on the orbit space V . li 6 = 0, then the result follows. Thus, we 
assume |^| > in C/ and show this implies that g has a free isometric local M x M action. 
Consider the 1-parameter family of metrics 

9s = T^*9V + + s-H)\ (2.4) 

for s > 0, with gi = g. Geometrically, this corresponds to rescaling the length of the fibers of vr 
and changing the horizontal distribution of tt, (when 0^0). Now it is a standard fact that the 
1-parameter family of 4-metrics 

gi=gs±u^dt^ (2.5) 

remains Ricci-flat for all s. This can be seen from standard formulas for Riemannian line bundles, 
c.f. [Bes, 9.36, 9G], [Kr, 16.1-3] or [A2, §1.2]. Thus, the metrics gs all satisfy the static vacuum 
equations 

urg = DgU 

from (0.1), with the same potential u. Equivalently, the conformal metrics gg = u^gs satisfy 

rs = 2{dloguf, (2.6) 

for all s, c.f. (1.4) The right side of (2.6) is of course independent of s. 

We claim that the metrics Qs on U arc all locally isometric. While this could be proved by a 
lengthy direct computation, we argue more conceptually as follows. Let e| be a local orthonormal 
frame for gg, determined as above for g. We then have ef = 61,63 = s~^63 while 63 varies in the 
plane < 62,63 > . The same relations hold w.r.t. gg- Recall also that the full curvature tensor is 
determined by the Ricci curvature in dimension 3. It then follows from these remarks and (2.6) that 
for each q e U and s > 0, there is a sectional curvature preserving isomorphism Fg : TqM TqM, 
i.e. 

kg{Fg{P))=Ki{P), 

where P is any 2-plane and Kg is the sectional curvature w.r.t. cjg. Clearly Fg varies smoothly with 
q and s. Using the expression (2.6), a result of Kulkarni [Ku] then implies that the metrics gg are 
locally isometric and hence so are the metrics gg. 

Let VL = dO he the curvature form of the line bundle vr. Then w.r.t. the metric g, \^\ = 
|^(ei, 62)1 = I < Vei62, 63 > |. The same equalities hold w.r.t. fig = d9g = s~^d9 and the gg metric. 
A short computation then gives 

\^s\ga as s ^ 00. (2.7) 

Hence consider the behavior of the metrics as s — > 00. We are then expanding or blowing up 
the metric in the fiber direction, at a given base point. Since the metrics gg are isometric, there are 
(local) diffcomorphisms 'ijjg such that 'ilj*gg converges to a limit metric g^o- At a given base point, 
the diffeomorphisms ipg expand or blow up smaller and smaller intervals of the parameter (p to unit 
size, giving rise to a limit parameter 4>oo- When ^ = 0, this is the only change; the limit metric ^oo 
is the same as g with the parameter ^ replaced by ^00 • (This is completely analogous to passing 
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from the flat metric on M? \ {0} to the flat metric on its universal cover, i.e. unwrapping the circles 
to lines). 

However, when 6^0, the 62 direction is also being expanded or blown-up in a similar way. The 
function u is left invariant under the family {ips}- 

It follows from (2.7) that the limit metric g^o is a static vacuum solution of the form 

5oo = 7rVoo + /^(#cx>)', (2.8) 

i.e. the 1-form 9^0 = in the limit. Further, since the 62 direction has been blown up, the function 
/oo varies only in the ei direction, i.e. /oo = foo{u)- 

Metrics of the form (2.8) are analysed in detail below. Referring to (2.10), let r = foo ■ u = h{u). 
In a possibly smaller open subset of U, we may invert h and write u = u{r), where r is a local 
coordinate on Voo- It is easy to see, (c.f. (2.12)-(2.13) below for example), that Qoc admits a non- 
vanishing Killing field d/dz, tangent to Foo but orthogonal to d/dr, and hence g^a admits a free 
isometric local M x M action. The metrics gg are all locally isometric and so the metric g = gi also 
has a free isometric local R x R action. 

■ 

Since the proof above is completely local. Proposition 2.2 holds locally, (in suitably modified 
form), even if (M, g) admits only a local or partial M -action. 

Static vacuum solutions admitting a free isometric local R X R action are completely classified; 
they are either flat or belong to the family of Kasner metrics, c.f. Example 2.11 below or [EK, 
Thm.2-3.12]. Such solutions do have Killing fields K which are not hypersurface orthogonal, i.e. 
dO in (2.3). For example, if d/dil^ and d/dz are standard generators of the (local) R x R action, 
then linear combinations such as K = 8/ dip + d/dz are non- hypersurface orthogonal Killing fields. 

Nevertheless, all such solutions do admit, of course, hypersurface orthogonal Killing fields and 
so may be written in the form (2.1). For the remainder of the paper, we thus assume that a Weyl 
solution has the form (2.1) locally. In addition, we will always work with the Z-quotient of the 
metric (2.1) and so consider Weyl solutions as warped products of the form V it will not be 

assumed in general that V is simply connected. 

Duality. Observe that Weyl solutions naturally come in 'dual' pairs. Namely the Ricci-flat 
4-manifold {N,gj^) has the form 

gN = gv + fd<t? + u^df, (2.9) 

and so both M„ = V Xf and Mj = V are static vacuum solutions on the 3-manifolds, 

with potentials u, resp. /. Consider the product of the lengths of circles, or equivalently, the area 
of the torus flber in N, 

r = f -u. 

This is a globally defined positive harmonic function on {V,gv)- To see this, on M 
we have 

= Au = Avu+ < Vlogf, Vtt >, 

and the same formula, with u and / reversed, holds on Mf. Hence 

Avfu = fAyu + uAvf + 2 < Vn, V/ >= 0. 

Charts. We now describe a collection of preferred charts in which to express the Weyl solution 
(M, g) ; this description is due to Weyl [Wl] . The surface V may be partitioned into a collection of 
maximal domains Vi on which the harmonic conjugate z of r is single- valued, so that F = r + iz is 
a well-defined holomorphic function from Vi into the right half-plane C"*" = {{r,z): r > 0, z G M}. 
One might also pass to a suitable cover, for instance the universal cover, of V to obtain a globally 
defined conjugate harmonic function, but it is preferable not to do so. 



(2.10) 
= M„, by (2.5), 
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Now each Vi may be further partitioned into a collection of domains Vij on which F is a conformal 
embedding into C"*", so that g\y.. = F*{dr'^ + dz^). We will thus simply view Vij as a domain in 
with gy a metric pointwise conformal to the flat metric dr^ + dz^. 

It follows that the corresponding domain Mij = Vij Xf is embedded as a domain = fi^ 
in endowed with cylindrical coordinates (r,z, (?!)). (6 G [0, 27r) with the background (unphysical) 
complete flat metric dr^ + dz'^ + r'^dcfP'. We note that all the data above are canonically determined 
by the two Killing fields on N and thus the coordinates (r, z, (j)) are called canonical cylindrical or 
Weyl coordinates for {M,g). Of course $7 is axially symmetric, i.e. symmetric w.r.t. rotation about 
the z-axis. 

To express the metric g\ci in these coordinates, the field equations (0.1) imply that the function 

u = logu (2.11) 

is an axially symmetric (independent of 0) harmonic function on C M'^; this again follows in a 
straightforward way from computation of the Laplacian of u and / on Mf and M„ as above. A 
computation of the conformal factor for the metric gv, c.f. [Wd,Ch.7.1], then leads to the expression 
of g in these coordinates: 

g = -u-2(e2^(dr2 + dz^) + r'^dcf^), (2.12) 

where A is determined by as a solution to the integrability equations 

Ar = r{u^ — u"^), Xz = 2ri'ri'z- (2-13) 

The equations (2.13) mean that the 1-form oj = r(f^ — v'^)dr + 2ri'rVzdz is closed on 

Conversely, given any axially symmetric harmonic function u on & connected open set Jl in M^, 
if the closed 1-form lo is exact, (for example if 7ri(Q n C"*") = 0), the equations (2.13) determine 
A up to a constant and the metric (2.12) gives a solution to the static vacuum equations with 
symmetry. (The addition of constants to or A changes the metric at most by diffeomorphism or 
homothety) . 

It is remarkable that solutions to the non-linear vacuum equations (0.1) can be generated in this 
way by solutions to the linear Laplace equation on M^. 

Remark 2.3. (i). The Levi-Civita solutions in §1 are all Weyl solutions. However, the expressions 
for the A1-A3 metrics in (0.4), (1.7), (1.8) and the Bl metric in (0.9) are not in Weyl canonical 
coordinates. Note that the quotients of the A2 metric discussed above are no longer Weyl solutions, 
although they could be considered as local Weyl solutions; the M-action is only locally defined on 
the quotients. 

(ii). There seem to be no known Weyl solutions which can not be expressed globally in the form 
(2.12). 

Fixed Point Set. The behavior of solutions at the part of dM where either one of the two 
or M actions on N has fixed points requires special considerations. This is of course the locus where 
n = or / = 0, and hence includes the part CtP[ A oi the z-axis A = {r = 0} in any canonical 
coordinate chart. It is not necessarily the case however that this locus is contained in A, c.f. the 
end of Example 2.10. 

Note that given any Weyl solution (2.12), any covering of \ A induces another solution of the 
form (2.12), but with ^ parametrizing a circle of length 27rfc. For the universal cover {k = oo), the 
(/(-circle is replaced by a line. In fact, (2.12) is well-defined when (j) runs over any parameter interval 
[0, 27ra). Observe however that any asymptotically flat Weyl solution must have a = 1, since the 
metric must be smooth near infinity. Thus, we will assume a = 1 in the following, unless stated 
otherwise. 
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Now suppose there is an open interval J in ^ such that the functions u and A, and hence the form 
(2.12) extend continuously to J. The form g then represents a continuous metric in a neighborhood 
of J if and only if the elementary flatness condition 



is satisfied on J. On intervals where (2.14) does not hold, the metric g has cone singularities, so 
that it is not locally Euclidean. Prom (2.13), it is clear that if A has a extension to J, then A 
is constant on J. However, such constants may vary over differing components of Cin A. This will 
be analysed further in Remark 2.8. 

For the remainder of §2, we assume that 

M = n, 

so that the Weyl solution is given globally in the form (2.12). Let I be the set where v = — oo, i.e. 
the Gs set in C given by 



where n runs over negative integers. It is usually assumed in physics that / is non-empty, although 
this need not be the case; this will also be discussed further below. The set / corresponds to the 
horizon S of the Weyl solution (M,g), since u = on I. This correspondence is formal however, 
since the geometry and topology of / C is very different than that of S C {M,g), c.f. most of 
the examples below. For the same reasons, although M = CI topologically, the metric boundary 
dM of {M,g) is (most always) very different than the Euclidean boundary of Q C M^. 

In the following, we discuss some of the most significant possible behaviors for the potential 
function u, and the associated Weyl solution, in order to illustrate the breadth of these solutions. 
The discussion is by no means complete or exhaustive. 

(I). Q compact. 

Let M = J7 be any bounded, C°° smooth axisymmetric domain (i.e. connected open set with 
smooth compact boundary) in and let 4> be any C*^'" function on dQ, > 0, a € (0,1), which is 
axially symmetric about the 2;-axis. For simplicity, assume that Q, fl C"*" is simply connected. Let 
v be the solution to the Dirichlet problem 

Av = 0,i'\dn = (t>- 

Then u is also axi-symmetric about the z-axis, and hence u generates a Weyl solution as in (2.12). 

Suppose that for a given A; > 1 a G (0,1), </> as above is (7^'" on dO., but is nowhere C^"*"^ on 
do,. Then v extends to a C'^'" function on the Euclidean closure H and hence, from (2.13), the 
function A in (2.12) is also uniformly bounded. This means that the metric g is quasi-isometric to 
the flat metric on and hence the metric boundary of CI w.r.t. the Weyl metric g is the same as 
its Euclidean boundary. Since v is not C^+^ anywhere on dCl, this solution M = Q.is maximal, i.e. 
admits no larger static vacuum extension; smooth solutions of the static vacuum equations are 
analytic. Further v is bounded, so that n = e'' is bounded away from 0, and hence the solution 
{Cl,g) has no horizon. 

As noted in §0, the presence of the boundary dM is physically assumed due to the presence of 
matter or field sources. Thus, at least when k > 2, the vacuum solution {M,g) can be extended to 
a larger space-like domain {M',g') D {M,g) with non-zero stress-energy T in M' \ M. 

On the other hand, if A; = above, then the geometry of the metric boundary {dM,g) will in 
general be very different than the smooth geometry of dCl in M^. Further, one can of course consider 
non-smooth domains O C in this situation. These remarks indicate that the structure of the 
metric boundary dM seemingly can be quite arbitrary. 



A = 



(2.14) 




(2.15) 



n 
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For the remainder of this section, we assume that Cl is non-compact in M^. The same remarks as 
above hold for non-compact domains with smooth (surface) boundary. Thus for example M = fl 
might have infinitely many ends if is non-compact, showing that the assumption (i) in Theorem 
0.3 is necessary. For simplicity, we only consider the following situation from now on. 

(II). Suppose 

dimndn < 1, (2.16) 

where the boundary is in the topology of R'^ and dim-j-i is the Hausdorff dimension. Thus dVL is a 
closed set of capacity 0, c.f. [H, Thm.5.14], and so in particular is a polar set. Clearly Cl = M^. 
(A). (Positive Case). Suppose that u is locally bounded above, i.e. 

sup < oo,Vx G (2.17) 

Bx(l) 

It follows, c.f. [H, Thm.5.18] that u extends uniquely to a globally defined subharmonic function 
on M^. Hence, one may use the value distribution theory of subharmonic functions on to analyse 
the geometry of Weyl solutions. 

The Riesz representation theorem c.f. [H,Thm.3.9], implies that any subharmonic function v on 
may be represented semi-globally, i.e. on B{R) = Bq{R) C for any R < oc, as 

i^ix) = - [ ^^d/x? + Hx), (2.18) 

Jb(r) F - ?l 

where d/j,^ is a positive Radon measure on B(R) called the Riesz measure of v and /i is a harmonic 
function on B(R); both d^ and h are axi-symmetric if v is. (A Radon measure is a Borel measure 
which is finite on compact subsets). For the moment, we only consider the situation where there 
exists K < oo, independent of R, s.t. 

/ ^^dii^ < K. (2.19) 

Jb{r) \x - ^\ 

In this case, one obtains a global representation of v as 

v{x) = - I -^^—dus + h{x), (2.20) 

7m3 \X - ^\ 

where djj.^ is a positive measure and h a harmonic function on M^. (In (D) below, we briefly discuss 

the situation where (2.19) is not assumed). In particular, if u is uniformly bounded above, say 
sup 1^ = 0, then the Liouville theorem for harmonic functions implies that h = 0, and one has the 
expression 

''i^) = - I Y^\dti^- (2-21) 
Jm? \x - K\ 

Note that since v is harmonic on $7, 

I C suppdiJL C (2.22) 
but in many situations, as will be seen below, the first inclusion is strict. 
(A)(i). Pure harmonic potentials. 

Suppose that z/ is a smooth harmonic function defined on all of M^, so that v = h in (2.20). It 
is clear that in this case S = in the sense that / = in M'^. Since ly is axisymmetric, v may be 
viewed as an expansion in Legendre polynomials, i.e. 

z. = ^afci?'=-Pfe(|), 

A;>0 

where R^ = + z^. Por instance, R ■ Pi(^) = z, R'^ ■ P2(^) = 3z^ - R^. 
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While these solutions arc defined on all of M^, no such solution gives a complete Weyl metric g 
on M = R^, by Theorem 1.1. For instance, for y = z, the Weyl metric is 

g = e-^'{e-'\dr^ + dz^) + r'^dcf^'^)), u = e^ 

Any straight ray in the (r, z) half-plane has finite length in this metric, except a ray parallel to the 
negative z-axis. The horizon S occurs formally at {z = — oo}, of infinite ^'-distance to any point in 

(A)(ii). Newtonian potentials. 

Suppose that /i = in (2.20), so that ly is the Newtonian potential of an axisymmctric positive 
mass distribution dfi as in (2.21). This situation corresponds exactly to the Newtonian theory of 
gravity, (or equivalently the electrostatics of a positively charged distribution). While there is a 
vast classical literature on this subject, we will only consider the most interesting situation where 

suppdfi = I, (2.23) 

so that v approaches — cx) on a dense set in supp d/j,. 
The following Lemma characterizes this situation. 

Lemma 2.4. Let d/i be an axi- symmetric positive Radon measure on M^. Then 

suppdjjL = I <^ suppd/j, C A. (2.24) 

Proof: Suppose first that supp d/j, is not contained in A. Since djj. is axially symmetric, part of 
supp dfi, namely the part not contained in A, is then given by a union of circles about the z-axis. 
Suppose first that there is a circle C which is an isolated component of supp d^, so that dn\c is a 
multiple of Lebesgue measure on C. This case has been examined in [W1],[BW], and we refer there 
for details. In particular in this case the potential u is bounded below on and near C, and hence 
supp dfi 7^ /. If C is not isolated, then using (2.21), the same reasoning holds, since the measure 
dfi is then even less concentrated on the circles. 

On the other hand, if supp dfj, C A, then dfi is a positive Radon measure on A. Standard 
measure theory implies that the upper density of dn w.r.t. Lebesgue measure dA at a e A, i.e. 
limsup^^Q ^ii^silii^ is positive, for Lebesgue almost all a G suppd^. From the expression (2.21), it 
is clear that for any such a, v{x) —oo as x — > a. This gives the converse. ■ 

For the remainder of the discussion in (A), we assume (2.23) holds. From the theory of subhar- 
monic functions on R^, the set I given by (2.15) may be an arbitrary set in ^4 C R^, i.e. a polar 
set in A. Since countable unions of polar sets are polar, note that / is not necessarily closed in 
A C R^. (For example, let {zj} G A be an increasing sequence converging to a limit point z, with 
say dist{zi, Zi^i) = and let dfi = X^2~*^2:., where is the Dirac measure based at Zi. Then 
/ = {zi} and supp d/x = {zi U z}). 

Given any a; G R^, let mx{r) be the mass of the measure djj, in the ball -B-r(r), i.e. 

mx{r) = / dii. (2.25) 

This is a non- negative increasing function on R"*", for any given x and the limit 

m = limr-tooi^x{i~) > 0) (2.26) 

is the total mass of dfi. This agrees, up to a universal constant factor, with the (ADM or Komar) 
mass in general relativity, when the latter is defined, and with (0.7) for solutions with pseudo- 
compact boundary. Note that one may have m = +oo. 

Lemma 2.4 and a standard result from potential theory, c.f. [H, Thm.3.20], characterize the 
possible Riesz measures satisfying (2.23). 
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Lemma 2.5. A necessary and sufficient condition that a positive Radon measure dfi is the Riesz 
measure of an axi- symmetric subharmonic function v on with sup u = and supp d^ = I is 
that supp djjL C A, and, for any given x & A, 

rrixir) 



i: 



-dr < oo. (2.27) 



It is easy to see that a Weyl solution (M, g) generated by a potential u as in (2.21) for which supp 
dn = I IS a compact subset of the axis A, is asymptotically flat, in the sense preceding Theorem 0.1. 
Further, the simplest or most natural surfaces enclosing any finite number of compact components 
of /, and intersecting A outside /, are 2-spheres in M. Of course if suppdjj, C ^ is non-compact, 
then the solution cannot be asymptotically fiat. A simple example is the solution generated by the 
measure 

where C, parametrizes A and dA is Lebesgue measure on A. Observe also that such solutions do 
not have pseudo-compact boundary. 

It is worthwhile to discuss some standard examples of Weyl solutions and their corresponding 
measures. 

Example 2.6. (i).(Curzon Solution). From the point of view of the Riesz measure, perhaps 
the simplest example is the measure d^ given by a multiple of the Dirac measure at some point on 
A, so that u = —m/R, R{x) = \x\, is a multiple of the Green's function on R^. This gives rise to 
the Curzon (or monopole) solution, c.f. [Kr,(18.4)], 
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= e2"^/-^[e-"^'^'/-^'(dr2 + dz^) + r^d(t>% (2.28) 



with u = e~"*/^. Here $7 = \ {0}, = {0}, and it is often stated that gc has a point-like 
singularity (monopole) at the origin. However the geometry of dM is very different than that of a 
point. Namely the circles about the z-axis have length diverging to infinity as i? — > 0. Thus, small 
spheres R = e about {0} become very long in the direction, and very short in the transverse 6 
direction, forming a very long, thin cigar. In particular, as a metric space, dM = M. This is the 
first example where dM is non-compact but pseudo-compact. Of course dM = S, so that {M,g) 
is complete away from S. A more detailed analysis of the Curzon singularity is given in [SS]. 

Note that one could not have solutions with both directions expanding at dM, so that area 
dM = oo, with dM pseudo-compact. This can be seen by use of minimal surface arguments, c.f. 
[G]. 

(ii). (Schwarzschild Solution). The Schwarzschild metric (0.4) is a Weyl metric, with measure 
dj^ = \dA on [— m, m], where dA is the standard Lebesgue measure on A. The resulting potential 
V in. (2.21) is the Newtonian potential of a rod on the z-axis with mass density |, given by 

-5 = \^-9{ ^R\Xfl^Z ^^ ^here i?| = + (. ± mf . (2.29) 

As mentioned before, the horizon E here is a smooth totally geodesic 2-sphere of radius 2m and 
dM = T,. 

The Weyl solution generated by the potential a - us, for us as in (2.29) with a > and a ^ 1, is 
not isometric or homothetic to the Schwarzschild metric. The associated Weyl metric is no longer 
smooth up to the horizon; in fact S is not even a 2-sphere unless a = 1. 
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Remark 2.7. More generally, consider any Weyl solution generated by a Riesz measure d/i satis- 
fying (2.23). Observe that / = ^, the length of the circles in the Weyl metric, stays bounded 
away from and oo on approach to supp d//, if and only if 

logr -C <v< logr + C, (2.30) 

for some C < oo, since v = logu. From the expression (2.21), this occurs only for the Schwarzschild 
potential us- Briefly, the reason for this is as follows. The estimate (2.30) implies that the potential 
V = v—vs is bounded and given by convolution of dist"^ with a signed Radon measure d\. However, 
as in the proof of Lemma 2.4, if v is bounded then one sees that necessarily d\ « dA and further 
the Radon-Nikodym derivative dX/dA is 0, a.e. (dA). In other words, any point of non-zero density 
for dX w.r.t. dA gives rise to approximating points on which u is unbounded. It follows that dX = 
0, and hence v = 0. 

In all other cases, where / ^ on approach to a region in supp dfi, the length of the circles 
goes to 0, and hence this portion of S is singular, of dimension < 1, while where / oo, the length 
of the circles goes to oo and the corresponding portion of S is singular and non-compact, (as in the 
Curzon solution). Note that if supp dfi is compact, then in all cases, dM is pseudo-compact. 

Thus among the Weyl solutions given by a Newtonian potential, only the Schwarzschild metric 
is smooth up to S. This gives a very simple (local) confirmation of Theorem 0.1 in this special case. 

Example 2.6.(iii). (Superposition/Multiple Holes). Subharmonic functions of the form 
(2.21) form a convex cone. In particular, one thus has a natural linear superposition principle for 
Weyl solutions. This feature is another remarkable property of Weyl solutions. 

For example, one may choose the measure d/Lt = ^dA on two, or any number of disjoint intervals 

{Ij} on the axis A, provided (2.27) holds. These correspond to solutions with 'multiple black 
holes', each interval Ij giving a component of S which is a 2-sphere of radius equal to the length 
of Ij. Although such solutions are essentially smooth up to S, they do not define smooth vacuum 
solutions on \ UBj. There are cone singularities, (called struts or rods in the physics literature), 
along geodesies (corresponding to A \ {Ij}) joining the 2-spheres of S, so that the metric g is not 
locally Euclidean along such curves. Thus, the elementary flatness condition (2.14) is violated on 
A \ Lilj. Nevertheless, the curvature of such metrics is uniformly bounded everywhere. These cone 
singularities must be considered part of dM, so dM in this case is a union of 2-spheres joined by a 
collection of curves and thus connected. 

Of course, the black hole uniqueness theorem. Theorem 0.1, also implies that such solutions 
cannot be smooth everywhere, when the number of intervals is finite. However, the proof of this 
result strongly uses the asymptotically flat assumption. It seems to be unknown whether there are 
any Schwarzschild type metrics with infinitely many black holes, i.e. metrics complete away from 
S and smooth up to E with S consisting of infinitely many 2-spheres, and which satisfy (2.17) or 
(2.21). It is natural to conjecture that such solutions do not exist, c.f. however the end of Remark 
2.8 below. 

Remark 2.8. Whenever supp d/j, C A is not connected but compact, there will exist such cone 
singularities on ^\ supp dfi. When supp dfi, or a sufficiently small smoothing of supp dfi, is 
interpreted to represent a matter source, then this statement corresponds exactly to the fact that 
there are no equilibrium (i.e. time-independent) many-body solutions in Newtonian gravity of this 
character, c.f. [Bn]. Even when supp dfi is non-compact and disconnected, this seems very likely 
to be true, c.f. the expressions for sums of Schwarzschild rods (2.29) in [IK, p. 336-337], which 
generalize to infinitely many rods. 

The components {Ci} of A\ supp dfi represent idealized matter sources (struts or rods) keeping 
the components of supp d/j, apart in static equilibrium. The cone angle a = is constant on each 
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C = Ci, and corresponds to a concentration of scalar curvature on C C M given by a multiple of 
the Lebesgue measure on C 

s = (1 — a)dAc, 

Thus the vacuum equations (0.1) are not satisfied across {Ci}. If a very small tubular neighborhood 
of radius r of such a rod is smoothly filled in with a perfect fluid source of constant pressure p and 
energy density fi, then one has the relation lirrir^orp = —limr^orp > 0. The effective mass of 
such rods is zero, i.e. they do not contribute to the gravitational potential v, c.f. [12] for a detailed 
discussion. 

By passing to covering spaces, it is always possible to create such cone singularities in Weyl 
solutions, even if none existed to begin with. For instance, for the Schwarzschild solution (0.4), with 
potential (2.29), take any covering, including the universal covering, ofE^\^ = (S'^\{aU— a}) x 
where {a, —a} are two antipodal points on S'^. This gives a solution whose metric completion has 
cone singularities along two (radial) geodesies starting at the antipodal points on = T, and going 
to infinity. 

Note that this discussion assumes that u is given by a Newtonian potential (2.21). In fact, there 
are Weyl solutions (M, g, v.) everywhere smooth up to the axis A, with S = dM disconnected, with 
no cone singularities or struts keeping the components of S apart. Namely, the Bl solution (0.9), 
dual to the Schwarzschild solution, has this property. 

Another, more remarkable, example is given in [KN]. These authors construct a Weyl solution of 
the form (2.12), which is complete away from S and smooth up to S, with S consisting of infinitely 
many Schwarzschild-like 2-spheres. In fact, the solution is periodic in the z-direction. This metric 
is not of the form (2.21), but is a limit of a sequence of solutions of the form (2.20), c.f. (D) below. 

Remark 2.9. If {M,g,u) is a Weyl solution of the form (2.12) with // = = logu. then the dual 
solution {M',g',f ), discussed in (2.9), is also a Weyl solution of the form (2.12), with potential i// 
= logf given by 

Vf = logr - Vu- (2-31) 

Hence if one potential is Newtonian, the dual one is not. Note that the sets lu^ If where and 

Vf are — oo are disjoint, with luL) If = A. Hence, if is a Newtonian potential with supp dp 
compact, so that the associated Weyl solution is asymptotically flat, then the dual Weyl solution 
is asymptotically cylindrical. 

Another example of a potential where both terms in (2.20) are non-trivial is the situation con- 
sidered (locally) in [GH] , where /i is a smooth axi-symmetric harmonic function defined on a neigh- 
borhood of supp dp, c.f. Remark 1.5. As vacuum solutions, these metrics cannot be complete away 
from S, as in the discussion on pure harmonic potentials. 

This completes our discussion of the Positive Case. 

(B) . (Negative Case). Under the assumption (2.16), suppose now instead that v is locally 
bounded below in R^, i.e. 

inf u > -oo,Vx G dn. (2.32) 

Then v extends uniquely to a globally defined superharmonic function on M^. Exactly the same 
discussion as in (A) above holds here, under the substitution v —v. (This corresponds to the 
transformation u following (1.5)). In this case, the Ricsz measure is a negative measure, 

so that one has solutions with negative mass. Note that here the potentials u oi u are unbounded 
above within supp dp, i.e. u or go to -|-oo. 

(C) . (Mixed Case). Next, one may superimpose Weyl solutions with positive and negative 
measures dp, i.e. consider u of the form (2.21), with dp a signed Radon measure. For example. 
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one may form dipole-type solutions with potential of the form u = + where and v- are 
(for instance) Curzon or Schwarzschild solutions of positive and negative mass placed at different 
regions on the axis. This gives for instance examples of asymptotically flat solutions with mass m 
assuming any value in M. 

More generally, since positivity is no longer assumed, the measure d/j, may be replaced by distri- 
butions, for example weak derivatives of measures. 

Example 2.10. (Multipole Solutions). As a typical example, one may take potentials corre- 
sponding to derivatives of the Dirac measure based at a point a G A, i.e. the multipole potentials, 

where Pn is the n*'* Legendre polynomial, or arbitrary linear combinations of such; c.f. [MF, 
p.l276ff]. 

Such potentials are limits of combinations of Newtonian potentials with positive and negative 
mass. Thus, it is reasonable to expect that there arc (Newtonian) equilibrium solutions, i.e. so- 
lutions with no cone singularities on the axis. This is proved in [Sz], where explicit equilibrium 
conditions are given. Note that one may have infinitely many multipole 'particles' in equilibrium. 

Another example is the potential of a dipole ring 



K^) = - / r^dC 
Jc F - 



z{x)_ 



where x = {r, z,(f>) and is the Lebesgue measure on the unit circle C = {r = 1} in the z = 
plane. Here z>(x) — > — cxd, as x — > C along the rays r = 1, z > 0. Hence, in this case, the set I = C 
is not contained in the axis A. 

(D). (Limits). Finally, one may consider potentials which are limits of potentials of the type 
(A)-(C) above, (besides those in Remark 2.9, Example 2.10). We consider just one important 
instance of this here. 

Example 2.11. (Kasner Metric). It is easily seen that the potential u = logr generates the flat 
metric 

g = dr^ + dz^ + dcj)^, 

on (M^)+, with potential u = r. Observe that since the 0-lines have constant length, the function 
r is now an affine (in fact linear) function on (M^)+. 

On the other hand, an equally simple computation shows that the potential v = a - logr, for any 
o G M, generates the metric 

g = r2"'-2"(dr2 + dz'^) + r'^'^" dcj)"^ , u = r". 

Equivalently, setting s = 

g = ds'^ + s'^dz'^ + s'^d(t>'^, u = s^, (2.33) 

where a = {2a-2)/{a-l+a-^), P = {2a~^-2)/{a-l+a'^),j = (a-l+Q-^)-^. Here s G M+,z G M 
and (f) G [0, 27r] or any other interval, including M, (by passing to covering and quotient spaces). 
These metrics arc all non-homothctic, provided a G [— 1, 0) U (0, 1]; a = gives the flat metric with 
u = 1 while a = — 1 gives the A3 metric (1.8). 

The potential v = a - logr can be considered as the limit 

u = limm^ooo\^s{''n^) — log2m\, (2.34) 

where usim) is the Schwarzschild potential (2.29) of mass m. Thus it is a limit of potentials of the 
form (2.18), where both terms are non-zero, (the harmonic term h is of course constant here). 
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These metrics are dual, in the sense discussed in (2.9) to the Kasncr (or Bianchi I) vacuum 
cosmological models, with homogeneous (flat) but anisotropic space-like hypersurfaces, c.f. [Wd, 
Ch.7.2]. It is easy to see that the Kasner metrics are the only Weyl solutions {M,g) which have an 
isometric R x M action, even locally. (The axisymmetric potential on $7 C must be invariant 
under an orthogonal M-action, hence giving a rotationally invariant harmonic function on R^. Thus 
the potential must be a multiple of log r). 

Consider these metrics on the quotient M = R+ x x S^.ln case a > 0, we have a < 0, ^ > 0, 
7 > and so S = dM. As in the discussion with the Curzon metric, the z-circles have unbounded 
length as s 0, so that dAI = R and the levels t = e, {t{x) = distg(x,dM)), are long, thin cigars. 
Thus, DM is non-compact, but pseudo-compact. The end of M \ J7 is obviously small. 

If a < 0, then a>0, /3>0, 7<0 and so S = 0, (it occurs at infinity), with dM = {pt}. In 
this case, the end of M \ ?7 is not small; the area growth of geodesic spheres is 0{r^~'^). 

However, none of these solutions are asymptotically flat even in a weak sense, except of course 
when a = 1. Namely, the curvature decays only quadratically in the ^'-distance to dM, i.e. 

\r\=0{t-^), 

and not any faster. Hence, the case a > shows that the conclusion of Theorem 0.3, (asymptotically 
fiat or small ends), cannot be strengthened to only asymptotically fiat ends, while the case a < 
shows that the assumption (ii) on u in Theorem 0.3 is necessary. 

Another metric of this (limit) type is that constructed in [KN], referred to in Example 2.8. This 
metric has the same asymptotics as the Kasner metric. 

3. Characterization of Asymptotically Flat Solutions. 

In this section, we prove Theorem 0.3. The proof of the first statement on finiteness of the 
number of ends is quite easy, so we begin with this. 

Throughout this section, let (M, g, u) be a static vacuum solution with dM pseudo-compact. We 
recall from §0 that M is connected and oriented. As in §0, let 

t{x) = distg{x,dM), 

and suppose U = t~^{0, Sg), so that dU fl M is compact. For r,s > So, let S{s) = t~^(s), A{r, s) = 
t~^{r, s) be the geodesic spheres and annuli about dU. It is important to note that neither S{s) nor 
A{r, s) are necessarily connected, even if M has only one end. (Of course if is a given end of M, 
then Se{s) = S{s) fl E must be connected for some sequence s = Sj —> oo). Let Sc{s) and Ac{r, s) 
denote any component of S{s) resp. A{r,s), so that S{s) = LiSc{s), A{r, s) = \jAc{r,s). Of course 
i is a proper exhaustion function on M \ U, so that these sets have compact closure in M. 

Let diam^Ac{r, s) denote the intrinsic diameter of Ac{r,s), i.e. the diameter of the connected 
metric space {Ac{r, s), g). 

Lemma 3.1. There exists a constant do < oo, independent of s, such that the number of compo- 
nents of A{^s,2s) is at most dg and 

diam^Ac{'^s, 2s) < do ■ s. (3-1) 

In particular, the manifold (M \ U,g) has a finite number of ends {-Ei}. 

Proof: Consider first the 4-manifold {N,g^), as in (0.2) which is smooth and Ricci-flat outside 
U = 7r^^(J7), where vr : ^ M is projection on the first factor. Since dU fl is compact, it 
follows from results of [Lu] that Lemma 3.1 holds on N, so that in particular N has a finite number 
of ends. Since N = MxS^,M\U also has a finite number of ends. 
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The choice of the time parameter on N defines a totahy geodesic embedding M G N and we 
have InIm = t where tN{x) = distN{x,dM),dM C —N. A geodesic ball or annulus in M embeds 
in the geodesic ball or annulus of the same size in A^. Hence (3.1) also holds for M. 

m 

Lemma 3.1 of course proves the first statement of Theorem 0.3. Observe that the estimate (3.1) 
is invariant under rcscaling of the metric g. 

For the remainder of the proof, we (usually) work with a given end E from the finite collection 
{Ei}. The main statement of Theorem 0.3 is that if 

^ -ds < oo, (3.2) 



areaSE{s) 

then the end E is asymptotically flat. The proof of this result is rather long, so we outline here 

the overall strategy. The asymptotic behavior of (E, g, u) is studied in general by examining the 
structure of the possible tangent cones at infinity, defined below. Basically, tangent cones at infinity 
fall into two classes, according to whether the asymptotic geometry near a given divergent sequence 
of base points is non-collapsing or collapsing, c.f. Lemmas 1.3-1.4. The main point is to prove that 
under the bound (3.2), all tangent cones at infinity are flat manifolds, and further that no collapse 
behavior is possible. Once this is established, the proof that {E, g, u) is asymptotically flat is 
relatively straightforward. 

Apriori, the end E may be very complicated topologically, for instance of infinite topological 
type; consider for instance that E might be of the form 5*00 x S^, where 5oo is any non-compact 
surface of infinite topological type and one end. A main idea is to use the behavior of the potential 
function u, in particular its value distribution theory, to control the topology and geometry of E 
in the large. We have already seen in §2 that the potential u controls quite strongly the geometry 
of Weyl solutions. Lemma 3.6 below is the key technical lemma which expresses this control for 
general static vacuum solutions (with pseudo-compact boundary). Further remarks on the strategy 
of proof precede the Lemmas below. 

The discussion to follow, until the end of Lemma 3.6, holds in general for ends E of static vacuum 
solutions with compact boundary. The estimate (3.2) will only be used after this. 

We now define the tangent cones at infinity of a given end E. (While this is a commonly used 
terminology, such limit metric spaces are not necessarily metric cones in general). 
First, we recall by Theorem 1.1 that there is a constant K < oo such that, Vx G M, 

\r\{x)<-^,\dlogu\{x)<-^. (3.3) 

The scale-invariant estimates (3.3) give quite strong initial control on the asymptotic geometry of 
{E,g,u) which allows one to get started. Observe that an immediate consequence of (3.1) and 
(3.3), by integration along paths in Ac{^s,3s), is the following Harnack inequality: 

= < K„ (3.4) 
mju 

where the sup and inf are taken over any component A^^s, 2s) and Ki is independent of A^^s, 2s). 

Let Xi be any divergent sequence of points in E, with ti = t{xi) oo. Consider the connected 
geodesic annuli Ai = Ai{K) = Ac{K~^ti, Kti),Xi G Ai, w.r.t. the rescaled or blow-down metric 

9i = ■ g; (3-5) 

here k is any fixed positive constant > 1. By the curvature estimate (3.3), the metrics gi have 
uniformly bounded curvature on Ai - the curvature bound depends only on K and k. Further, by 
(3.1), the diameter of A^ w.r.t. gi is also uniformly bounded. 
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Hence, if the sequence is non-collapsing, i.e. if there is a lower volume bound volgAi > I'o'tf, for 
some Vo > 0) (equivalent to volg^Ai > by scaling), then Lemma 1.3 implies that a subsequence 
of the pointed sequence {(^i, 5^, Xj)}, converges smoothly, away from the boundary, to a limiting 
smooth metric (^ooC/^); 5oo; a^oo)- The limit is a solution of the static vacuum equations; as noted in 
Lemma 1.3, the potential u is renormalized to Ui = u/u[xi), so that the limit potential Uoo satisfies 
'^oo{xoo) = 1; c.f. also (3.4). Choosing a sequence nj 00 and a suitable diagonal subsequence, 
gives the maximal static vacuum solution (^ooj 5oo) ''^ooj a^oo)- Observe here also that the estimate 
(3.1) implies that dA^o = {pt}. 

On the other hand, if the sequence {Ai,gi) is collapsing, in the sense that volgAi << tf, as 
ti 00, (equivalent to volg^Ai 0), then Lemma 1.4 implies that Ai is a Seifert fibered space or 
torus bundle over an interval. As discussed there, one may then pass to Z or Z Z covers Ai{K,) 
of Ai^K) (or more precisely smooth interior approximations to Ai{K)) to obtain a non-collapsing 
sequence {Ai{K), Qi, Xi), smoothly convergent to a limit {Aoo{k), g^o, x^o)', {gi here is lifted to the 
cover ylj, as is Xj). As above, one may then choose a sequence Hj — >■ 00 and pass to a diagonal 
subsequence to obtain a maximal limit (^00 5oo) '"oo)- This limit static vacuum solution has an 
isometric M action, or M0M action in the case of a rank 2 collapse. Hence by Proposition 2.2 it is a 
Weyl solution. In the latter case, the solution is then a (possibly flat) Kasner metric, c.f. Example 
2.11. Topologically, the limit Aqq is a trivial M, (or M © M), bundle over a surface F, (or interval), 
again by Proposition 2.2. 

As in §2, we will always work in a Z-quotient A^ of Aoo, (or Z © Z quotient in the case of 
rank 2 collapse), and finite covers Ai{K) converging to Aoo(/^)• For any fixed k > 0, the manifolds 
{Ai{n),gi,Xi) thus have uniformly bounded curvature and diameter, a uniform lower bound on 
their volume, and converge smoothly to the limit (^oo(«), ^oo) C (^00)500)- The limit has a free 
isometric or x action, and so in particular is an or bundle. Hence ^j(/«) is also 
topologically an or bundle, (although not metrically). 

To be definite, the finite covers are chosen so that the length of the factor or factors at the 
base point Xj G Ai{K) converge to 1 in the limit. 

Recall by Lemma 1.4 that the inclusion map of the fibers induces an injection on vri. The coverings 
Ai^n) arc obtained by taking large finite unwrappings of the or fibers, (corresponding to 
taking subgroups of 7ri(5^) or 7ri(T^) of large but finite index). All finite covering spaces of 
or are still or T^, and hence we may, and do, choose the unwrappings so that, as smooth 
manifolds, 

Ai{K) = Ai{K), 

for any k > 0. In the limit, the unwrapping of the collapse thus just corresponds to expanding the 
length of the collapsing factor (or factors), preserving the holonomy, if any, of the bundle; 
compare with the proof of Proposition 2.2. 

The limit spaces (^oo 5oo) 'Jioo) a;oo) or (^00 ^oo ^^ooi a^cxa) constructed above are called tangent 
cones at infinity of {E,g,u). Note that such tangent cones are only attached to some subsequence 
of a given divergent sequence of base points {.Xj}. Hence, apriori, the tangent cones at infinity could 
be highly non-unique as Ricmannian manifolds. In general, there may be no relation between the 
geometry of different tangent cones based on (subsequences of) distinct divergent sequences {xi}; 
for example, tangent cones based on sequences with t{xi) = 2' and t(xj) = 2' . The tangent cones 
only detect behavior of the end E in ^j-bounded distance to the base points Xj. 

On the other hand, since tangent cones at infinity attached to any divergent sequence always 
exist, for any s sufficiently large, say s > So, the geometry of (^c(|) '^),9s) or {Ac{^, 2),5s) is always 
close to that of some tangent cone at infinity ^00 or ^00 ■ Further, by construction, the tangent 
cones are always connected and, since M is oriented, so is each tangent cone. 

The following lemma is a typical application of the use of tangent cones at infinity. 
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Lemma 3.2. Suppose the curvature r decays faster than quadratically in the end {E,g,u), i.e. 

M(x) < (3.6) 

where £{t) as t ^ oo. Then there is a compact set K C E such that E\K is dijfeomorphic 
either to \ B or to (E? \ B) x S^, where B is a 3-ball, (resp. a 2-hall), i.e. E is of standard 
topological type. Further, the annuli Ae{\s,2s) are connected, for all s > Sq, for some Sq < oo. 

Proof: By the preceding discussion, the condition (3.6) is equivalent to the statement that ah 
tangent cones at infinity (^oo^oo) or (^oo^oo) of E are fiat, (as well as connected and oriented). 
The two possible conclusions of Lemma 3.2 correspond to the two possibilities of non-collapse and 
collapse in the formation of the tangent cones. 

Suppose first that g is non-collapsing on E, i.e. there exists Uo > such that volgAc{'^s,2s) > 
Vo ■ s^, for all s large, and all components Ac. Recall that dA^o = {pt} in this situation. It then 
follows that for s sufficiently large, each Ac{^s,2s) is diffeomorphic, and almost isometric to the 
standard flat annulus A = r~^{^s,2s) in M^, r{x) = \x\, (away from the boundary). In fact, each 
tangent cone at infinity Aqo is isometric to \ {0} in this situation. Here we are implicitly using 
the fact that the only complete oriented flat 3-manifold with an isolated singularity is \ {0}, 
c.f. [AC] for example. Similarly, a smooth approximation to Sc{s) is diffeomorphic and almost 
isometric to S'^{s) C M^. 

By the isotopy extension theorem, these diffeomorphisms from Ac{^s, 2s) to the standard annulus 
may then be assembled to a global diffeomorphism, and almost isometry, of E \ K into M.^ \ B, 
for some compact set K C E. We refer to [AC, Thm.1.18] for the proof of these statements, (in a 
slightly different but equivalent form), which are now quite standard. The main point is of course 
that since the family of annuli Ac{^s,2s) as s varies is topologically rigid, i.e. one has a unique 
topological type, there is no value of s at which the topology can change or bifurcate. 

Suppose on the other hand that g is collapsing on Ac{^Si,2si), for some sequence Sj oo, 
and some sequence of components Ac. As discussed above, one may then pass to suitable covers 
Ai = Ac{^Si,2si) so that, in a subsequence, {Ai,gi) is diffeomorphic and almost isometric to its 
limit Aryo{^, 2) C Aoo- The maximal limit Acx^ is a flat manifold with either a free isometric or 
X action. Hence there are two possibilities for A^o, namely either V x or M.'^ x x S^, 
where K is a flat 2-manifold and the metric is a product metric on each factor. In the former 
case, the diameter estimate (3.1) imphes, (as in the non-collapse case above), that F is a complete 
flat cone, possibly with an isolated singularity at {0}. Hence, although these two possibilities 
for the limiting metric of ^00(5; 2) are distinct, both are the same topologically, i.e. Aoo{^,2) is 
topologically I x x S^. 

Now recall from the discussion on tangent cones that A^ is diffeomorphic to Ac{^Si,2si); metri- 
cally Ai approximates one of the types of flat manifolds above, with factors shrinking to very 
short circles. In both cases, (a smoothing of) Sc{si) is diffeomorphic and almost isometric to a flat 
torus T^. 

In particular, the topological type of A^ is distinct from that of the annuli Ai above in the non- 
collapse case, which are topologically always of the form / x 5^ C M'^ \ {0}, (for any choice of base 
point and component). This implies first that the family {Ac{^s,2s)} must be collapsing for all s, 
as s — > 00, and all components Ac. Second, the topological type of the annuli Ag, and hence that 
of As = Ac{^s, 2s) is unique, and given for all s large and all chy I x x S^. Use of the isotopy 
extension theorem in the same way as above then proves that the end E itself is diffeomorphic to 

(M2 \5) X S^. 
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Remark 3.3. (i) It is easily seen from the vacuum equations (0.1) that the condition (3.6) is 
equivalent to 

\Vlogu\{x) < (3.7) 

as i — > oo, c.f. also the proof of Theorem 1.1. 

(ii) . In the context of Lemma 3.2, the tangent cones at infinity (^oo,5oo) may not be unique 

up to isometry in the collapse case, and so may vary within the moduli space A^i of flat product 
metrics of the form V x or within the moduli space M2 of flat product metrics of the form 
M+ X 5^ X S^. 

Note that the moduli space Mo of flat metrics on or \ {0} is just one point, (c.f. again 
[AC] for the latter statement for example). Similarly, by the normalization preceding Lemma 3.2 
that the factors have length 1, the moduli space A^2 -^2 normalized in this way is also just 
one point. The moduli space M'l C A^i where the factor has length 1 is naturally identified 
with R"^, parametrized by the cone angle at {0}. 

Observe however that these two moduli spaces A4'i and M2 are disjoint; they cannot be connected 
(or even approximated) by a curve of flat metrics. Now the geometry of the annuli {Ae{^s, 2s), gs) 
varies continously with s. By the remarks preceding Lemma 3.2, this induces a continuous variation 
of the possible tangent cones {Aoc,goc) in -^1 or A42- Hence, on a given end E, one cannot obtain 
two different tangent cones, one of the form V x and another of the form R+ x S"^ x S^; c.f. also 
the proof of Lemma 3.7 below. 

(iii) . We will need a slight generalization of Lemma 3.2 for the next lemma below. Thus let 7(5) 
be any properly embedded curve in E with ^(7(5)) — > cx) as s ^ 00, and suppose (3.6), (or (3.7)), 
holds in the balls B^(^g^(6 ■ t{'y{s))), for some fixed 5 > 0. Then the conclusion of Lemma 3.2 also 
holds. 

To see this, consider the blow-downs gg = t~'^{y{s)) ■ g, based at 7(5), and the associated tangent 
cones at infinity. The scale invariant condition (3.6), when applied to B^(^g-^{5 ■t{'y{s))), implies that 
all such tangent cones are flat in {Bx^{d), goo), and thus flat everywhere in their maximal domain 
Aoo, by the fact that smooth solutions of the vacuum equations are real analytic. The proof then 
proceeds exactly as in Lemma 3.2 

To prove that the estimates (3.6) and (3.7) do in fact hold on E, we need to understand in more 
detail the value distribution of the potential u. The main result needed for this is given in Lemma 
3.6, and then (3.6)-(3.7) follow rather easily in (3.17)-(3.18) below. However, some preliminary 
results are required for the proof of Lemma 3.6. The main difficulty is that u may not, in this 
generality, be a proper function onto its image, (c.f. the remark following the proof of Lemma 3.4). 
Lemma 3.4 below is a slightly weaker substitute for this property. 

Let U = t-^{0,So) be a neighborhood of dM as in the beginning of §3. Let 7(t) be a maximal 
flow line of Vn. We will say that 7 is divergent if 7 does not intersect U at two different times, i.e. 
if 7(r) exits U at some time, then 7(r) never reenters (a possibly distinct component of) U, and 
if further 7(t) does not terminate at a critical set of n in M \ [7 as r — > ±00. It follows that if 7 
is divergent, then 7 is complete in at least one direction, (r +00, r — >■ —00 or both), and in any 
such direction, j{t) diverges to infinity in M. Since the potential u has no local maxima or minima 
in M, the set of flow lines terminating on a critical set of u in M \ [7 is a closed set of measure 
in M. This follows for instance from the fact that the measure \du\dA, where dA is Lebesgue 
measure on the level sets of «, is preserved under the gradient flow of u, and this measure tends to 
on approach to critical points of u. Thus among the set of flow lines not joining points of U, the 
divergent flow lines are generic in terms of measure on M. Of course the flow lines are curves of 
steepest ascent for w as r increases, and of steepest descent for u as r decreases. 
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Lemma 3.4. There exists a compact set K C M , with U G K, such that any divergent flow line 
7(t) intersects K. 

Proof: Suppose that this were not the ease, so that there exist, necessarily complete, flow lines 
7(t),t G M, which do not intersect a given K D U. We may choose K sufficiently large so that 
7(r) is then contained in a fixed end E C M \ U, since there are only finitely many ends. 

Let Aj{t) = y4^(|t(7(T)), 2t(7(T))) be the component of the geodesic annulus containing the 
base point j{t) and let E"^ be the part of E swept out by such annuli, -B-y = Ut-Aj^t) C E. As 
r — > cx), the function u(7(r)) is monotone increasing. 

If n(7(T)) increases to +00, then n — > 00 uniformly as r — > 00 in E^ = UT->o^7(r) C E, by 
the Harnack estimate (3.4). Since E is an end, there exists some sequence tj 00 such that, 
for tj = i(7(Tj)), the spheres S^{tj) C A^{Tj) satisfy S^{tj) = SE{tj), i.e. the spheres SE{tj) 
are connected, and hence Aj{Tj) = AE{Tj). Thus, u becomes uniformly unbounded on AE{Tj), as 
j 00. However, as r ^ —00, the curve 7(r) also diverges to infinity in E and n(7(T)) is decreasing, 
(and so in particular bounded), as r ^ —00. This contradiction implies that u(7(t)) u'^ < +00, 
as r — > +00. Of course w+ > 0. 

Suppose first, (for simplicity), that limsupEU = i.e. limt-^oofn{i) = 'i^^; where m{t) = 
sups^(t)'f-''- The maximum principle for the harmonic function u implies that for t sufficiently large, 
the function m(t) is either monotone increasing or monotone decreasing in t and hence approaches 
the value u'^ast^ 00. Consider the annuli A^{t) in the scale gt = ■ g, t = t{'y{T)), as in (3.5). 
Any sequence r, — > cxd has a subsequence such that the corresponding annuli {A^{Ti),gti) converge 
to a limiting domain Aao{^, 2) in a tangent cone at infinity (^00, 5oo)) or (^00, fl'oo)) passing to covers 
as described above in the case of collapse. By construction, we then see that the potential function 
Uoo for this limit static vacuum solution achieves its maximal value n"*" > at an interior point. 
Since ■Uqo is harmonic, the maximum priniple implies that Uoo = u'^, and hence by the vacuum 
equations (0.1), the limit (^ocffoo) or (^oc^oo) is flat. This argument holds for any subsequence, 
and since the convergence to the limit is smooth, we see that 

■ '^U^(it,2t) ^ 0' ^ = ^(7(-^)), as 

by the scale-invariance of this expression. 

It follows from Lemma 3.2 and Remark 3.3(iii) that the end E is topologically standard, and 
the annuli AE{^t,2t), t = t(7(r)), are connected in E, for all r sufficiently large. Prom the prior 
argument, this implies in particular that u ^ uniformly at infinity in E. 

However, as before, as r — > —00, 'u(7(t)) is monotone decreasing to a value U- > 0. It follows 
that u'^ = U-. This is of course impossible, and shows that 7(r) must have exited E at some 
(negative) time. 

Thus it remains to prove that L = limsupEU = u^. Since the annuli AeItj) above are connected, 
the Harnack inequality (3.4), together with the fact that u'^ < 00, implies that L < 00. Now choose 
points Xj G S^{tj) such that u{xj) — > L. As above, the functions ^^^(rj) have a subsequence 
converging to a limit harmonic function u^q on a tangent cone at infinity based at 
Then as before u has an interior maximum at x^ and hence Uoo = L\ this gives L = . 

■ 

It follows from Lemma 3.4 and the discussion preceding it that any maximal flow line of 
intersects an apriori given large compact set K C M, except those exceptional flow lines which 
start or end at a critical point of u far out in M. In particular, a set of full measure in any given 
level set L of u may be connected to points in K by flow lines of Vu. In this sense, u is 'almost 
proper', in that it behaves almost like a proper function in terms of the gradient flow. 

Observe that this does not necessarily imply that the level sets of u are compact, i.e. that u is 
proper. For instance, the Weyl solution generated by the dipole potential z/ = z/+ + z/_ considered 



25 



in §2(IIC) satisfies (3.6), (it is even asymptotically flat), but the 0-level of v is non-compact if 
and z/_ are chosen so that the mass is 0. In this example, the only divergent flow lines of Vu are 
the two ends of the z-axis. 

Next, as in §2, let 

V = logu. 

The following result is quite standard. 
Lemma 3.5. On {N,g]s[) as in (0.2), with Riemannian metric, we have 

Ajv|Vz/| > 0. 

Proof: This standard estimate is a simple consequence of the Bochner-Lichnerowicz formula 

iA|Vi/|2 = |D^z^p+ < VAi^,Vz^ > +r(Vi/,Vi/), 

on {N,g]\[), where we have dropped the subscript N from the notation. Since Ai/ = 0, by (1.6), 
and since {N,gN) is Ricci-flat, this gives 



One computes 

and, by the Cauchy-Schwarz inequality |V|Vz^pp < 4|Z)^i/p|Vi^p, so the result follows. 



A|Vz/| = i|Vz^|-^A|Vi^|2 - i|Vi/|-3|V|Vi/|2|2, 



Lemmas 3.4 and 3.5 lead to the following key result relating the behavior of \Vu\ to the area 
growth of geodesic spheres. This result is a straightforward consequence of the divergence theorem 
for proper harmonic functions u on manifolds of non-negative Ricci curvature. Lemma 3.4 allows 
one to remove the assumption that u is proper. 

Lemma 3.6. There is a constant C < oo such that for any component Sc{s) of S{s) C M, s>l, 

sups^(^g^\Vu\ < C ■ areaSc{s)~^ . (3-8) 

Proof. We work on the Riemannian 4-manifold [N^ gi\f) until the end of the proof. Let Ac = 
Ac{s) = Tr"^{Ac{^s, 2s)) and Sc = Sc{s) = tt^^{Sc{s)), where tt : N ^ M is projection on the first 
factor, with fibers. From the coarca formula, we have 

/ |Vz/|^ = / / \S/i^\davdv, (3.9) 
J Ac Jv J L^nAc 

where Ly is the v-level set of in iV and the outer integral in (3.9) is over the range of values in 

R of in Ac- Now as remarked following Lemma 3.4, up to a set of measure in L^, D Ac, all 
points in the set (Ly H A^.) \ Zy may be joined by flow lines of Vz^ to points in a fixed bounded 
hypersurface T in K, independent of v, s; here K is the compact set from Lemma 3.4, and T = dK 
for instance. Hence, by the divergence theorem applied to the harmonic function v on iV, 

/ iVi^l < / |Vz/| < ci, (3.10) 

for some ci < oo, independent of s and Ac- 

Now is subharmonic on {N,gjsj) by Lemma 3.5, and diam^j^Ac < c- s, by (the proof of) 

Lemma 3.1. A standard sub-mean value inequality for manifolds of non-negative Ricci curvature, 
c.f. [SY,Thm.IL6.2], then gives 



supa rJVul'^ < I (3.11) 

' - vol A J Ac 
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Hence the estimates (3.9)-(3.11) imply 

supg^^^^\Vi'\'^ < C3 • osc^^u ■ (volAc)"^. (3-12) 

To estimate the right hand side of (3.12), again by (the proof of) Lemma 3.1, osc u < 04 - sup\\/u\s 
on Ac. Further, we claim that 

for some constant c independent of s and Ac- To see this, by scale-invariance, it suffices to prove that 
sup^ (1)1^^1 — (i)l^'^l w.r.t. the rescaled metrics Qg = s~'^g. By the curvature and diameter 

bounds (3.3) and (3.4) and Lemmas 1.3 and 1.4, the metrics {Ac{l),gs) form a compact family of 
metrics in the C°° topology, unwrapping in the case of collapse. Thus, one has uniform control 
on the metrics gs on Ac{l). Similarly, when normalized if necessary by additive and multiplicative 
constants so that sup^ ^-^-^ly = sup^ ^-^^[Vz^j = 1, the positive harmonic functions on (A(.{l),gs) 
also form a compact family of functions in the topology, i.e. a normal family. This follows by 
the Harnack estimate (3.4) and the Harnack principle (elliptic regularity) for harmonic functions, 
c.f. [GT, Thm 2.11, Ch. 8]. This compactness of the metrics and functions from elliptic theory 
proves the claim above. 

Similarly, the metric compactness above also implies there is a constant c < 00 such that 

■ areaSc{l) < volAc{l) < c ■ areaSc{l), 

w.r.t. the metrics gs- (This estimate can also be derived directly from the Bishop-Gromov volume 
comparison theorem). Rescaling back to the metric g then gives 

c^"*" • s ■ areaSc{s) < vol Ac < C5 • s • areaSds), (3.13) 

for some constant C5 < 00. Note that by definition, 

areaSc{s) = / udA, (3-14) 

JSc(s) 

and the same for volAc- Hence by (3.4), the estimate (3.13) holds also for Sc{s) and Ac in place of 
Sc{s) and Ac. 

Thus, by combining these estimates above, (3.12) gives 

supg (s)IV^'l < cg • areaSc{s)~^ . 
Using (3.14) and (3.4) again, this estimate implies (3.8). 

■ 

We are now in a position to begin the proof of Theorem 0.3 itself. Observe that the previous 
results in §3 have not used the assumption (3.2), nor the assumption (ii) in Theorem 0.3 that u 
does not approach everywhere at infinity in E. Only the assumption that dM is pseudo-compact 
has been used. Hence, at this stage, we do not even know that E has finite topological type. 
The main point initially is to prove that the estimates (3.6)-(3.7) above do hold on E under these 
assumptions. 

Recall that we have Se{s) = US'c(s), for Sc{s) C E. Each geodesic ray a{s) in E, i.e. an integral 
curve of Vt, with a{s) G Se{s), determines a component Sa{s) = Sc{s) s.t. a{s) G Sc{s); the union 
of such components sweep out a part E^^ of the end E. Of course E is the union of E^^ among all 
(non-homotopic) rays a. 

Prom Lemma 3.1 and from the obvious areaSE{s) = X] areaSc{s), we have 

/oo poo 
areaSE{s)~^ds < 00 4^ / areaSa{s)~^ds < 00, 
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for some geodesic ray a C E. (Here the integrals start at some fixed value s > So > 0). 
Hence, under the assumption (3.2), we have 

/oo 
areaSa{s)~^ds < K < oo, (3.15) 

for some ray a <Z E and constant K. 

By integrating along the curve cr, (3.15), Lemma 3.6 and the Harnack estimate (3.4) imply that 
u is uniformly bounded in E^. In fact, we claim that 

■Uoo = limt(^-)^^u{x) < oo (3.16) 

exists, where the limit is taken in E^r- To see this, let 7(5) be any "quasi-geodesic" in E^^, i.e. 7 is 
a smooth curve with 7(5) G S'o-(s) and Idj/dsl < Ci, for some Ci < 00. By (3.8) and (3.15), we 
then have J°° \du{'y{s))\ds < C ■ Ci ■ K < 00, and so ^(7(51)) — u{^{s2)) — > whenever si, S2 — *■ 00. 
Hence the limit itoo(7) is well-defined. The diameter estimate (3.1) implies that all points in E^- lie 
on quasi-geodesics, (with a fixed Ci), in E^^, starting on Sa{so)- Further, the limit 'Uc!o(7) is clearly 
independent of 7, since for instance (3.8) and (3.15) imply that 

on some sequence sj 00. Hence (3.16) follows. 

Next, since E is an end, there exists some sequence tj 00 such that the geodesic spheres 
SE{tj) are connected, and hence Sf^itj) = SEitj)- By (3.16), u\g^(^f.-^ — > Uqc as tj 00. The 
maximum principle applied to the harmonic function u thus implies that u\AE(tj,tk) ~^ ^00 whenever 
tj^tk cxD. Thus, we see that (3.16) holds where the limit is taken in the full end E, and not just 
in E„. 

Now we use the assumption (ii) of Theorem 0.3, which says that u{xj) > > 0, for some 
constant Uq > and some divergent sequence Xj G E. It follows from this and the existence of the 
limit (3.16) in E that 

Uoo > 0. 

Hence, wc may, and will, renormalize the potential function u of the static vacuum solution (M, g, u) 
so that, on E, 

limt{^)_^^u{x) = 1. (3.17) 

The estimate (3.17) essentially immediately implies the scale-invariant estimates 

supSe{s)\A « s^^P5b(s)|Vu| « as s ^ 00, (3.18) 

strenthening the bounds (3.3). For as discussed following (3.3), (3.18) is equivalent to the statement 
that all tangent cones at infinity of E are flat, (c.f. also the proof of Lemma 3.2). But, as noted 
above, all tangent cones at infinity are static vacuum solutions, and (3.17) implies that the limit 
potential u^o satisfies u^o = 1- Hence, the static vacuum equations (0.1) of course imply the limit 
metrics are flat. 

Lemma 3.2 now determines the topology of the end E, as one of two (standard) alternatives, 
according to non-collapse or collapse behavior at infinity. 

Before proceeding to the analysis of these cases, note that (3.18) implies that the metrics g and 
g = u^ ■ g from (1.4) are quasi-isometric on E, and almost isometric near infinity. Since g has 
non-negative Ricci curvature on E, standard volume comparison theory implies that the area and 
volume ratios 

areag{S{s)) volg{B{s)) 
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are monotone non- increasing in s. Hence, their hmits at s = oo exist, and by the equahty of g and 
g at infinity, the limits at s = oo of the ratios in (3.19) w.r.t. the g metric and ^'-geodesic spheres 
and balls also exist, and equal the g limits. 

To proceed further, we now separate the discussion into non-collapse and collapse cases. 

Case A. (Non-Collapse). 

Suppose that E is non-collapsing at infinity, i.e. by the remarks above, 

limsups^oo — ^ = lims^oo — ^ > 0, (3.20) 

where v{s) denotes volg{BE{s)). This implies, via (3.13), that areaS{s)~^ < c- s~^, and hence, by 
Lemma 3.6, 

sups(^g-j\Vu\ < cs~'^ . (3.21) 

The volume condition (3.20) implies that all tangent cones at infinity A^o of {E,g) exist, without 
passing to covering spaces, and, as discussed above, are fiat solutions of the static vacuum equations 
with potential Uqo = 1. As discussed in the proof of Lemma 3.2, the tangent cone at infinity is 
unique, (up to isometry), and given by \ {0}, and further £^ \ is diffeomorphic to \ B, for 
some compact set K C E. The blow-down metrics gg = • g on all annuli Ae{^s, 2s) converge 
smoothly to the flat metric on ^(1,2) C M^, uniformly as s — ^ oo, and hence there are local 
(harmonic) coordinates on ^£;(^,s, 2s) in which g has the expansion gij = ^ij+^ij, where |7jy(-T)| — > 
uniformly as s — > oo. Again as discussed in the proof of Lemma 3.2, these local coordinates, e.g. 
on {Ae{2^~^, 2*+^)}, i > 0, may be assembled into a global chart, mapping E\K onto \ B, c.f. 
[AC] for further details if desired. With respect to such a chart, the metric gij has the form 

gij = Sij + -fij, (3.22) 

on all of \ i^, with \^ij\ uniformly at infinity in E. In other words, g is C° asymptotic to 
the flat metric at infinity. 

To prove that the metric on E is (strongly) asymptotically flat, as defined preceding Theorem 

0. 1, consider again the metric g = ■ g. Prom (1.4) and (3.21), the curvature of g decays as 

\f \ < cr^, 

as t — > oo. (Note also that t and t are approximately equal for t large). 
It follows that the expansion (3.22) may be improved, for g, to 

gij = Sij + 0{t-^), (3.23) 

in a suitable (harmonic) coordinate chart. We refer to [BKN] or [BM] for instance for further 
details here. Briefly, elliptic regularity theory applied to the equations (1.4)-(1.5), together with 
the curvature decay above, implies that the 2"'' derivatives of the metric g in the coordinate chart 
decay as t~^, so that the metric g decays to the fiat metric at a rate of t~^. Hence 

gij = u-%j = (1 + 2v)5ij + 0(i-2), (3.24) 

where v = 1 — u. Here we are using that fact that since \Vu\ = 0{t~^), u has an expansion of the 
form u = 1 +0{t~^). Further, since log u is harmonic w.r.t. g, the decay (3.21) and (3.23) implies 
that Aflogu = 0(t~^) for t large, where Aj is the fiat Laplacian on M.^. This means that u has an 
expansion w = 1 — y + 0(t^^), where m is the mass of E defined in (0.7). 

In particular, these estimates show that the end E is asymptotically fiat in the sense preceding 
Theorem 0.1. 

Note that, to first order in t~^, the function v = l — u corresponds to the Green's function in M.^, 

1. e. the fundamental solution of the Laplacian, weighted by the mass m. It is of course possible to 
have m = 0, as for instance for the dipole-type Weyl solutions in §2(IIC), or also m < 0. Further, 
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since u has been normalized so that u ^ 1 at infinity in E, the expression (0.7) for the mass is 
equivalent to the usual definition 

mE = ^ I <Vu,Vt> dA, (3.25) 

4^ JsEis) 

where s is sufficiently large so that Se{s) n dE = 0. This is because the expression (3.25) is 
independent of s, since u is harmonic, and the fact that it is asymptotic to the expression (0.7) as 
s — > oo. 

This completes the analysis of Case A. 
Case B. (Collapse). 

Under the standing assumption (3.15), suppose that the end E is collapsing at infinity, i.e. 

limsups^oo-^ = lims-,oo^^^ = 0. (3.26) 

We will prove that this situation is impossible. The results preceding Case A remain valid, so 

that (3.17)-(3.18) hold, all tangent cones at infinity Aq^ of E are flat products of the form V x or 
M"*" X 5*^ X 5^, where V is a flat 2-dimensional cone. Further £^\i^ is diffeomorphic to (M^ \iJ) x S^, 
for some compact set K C E. 

By (3.15) and Lemma 3.6, we have 



sups{s)Nu\{s)ds <Ki<oo. (3.27) 



The main point is now to show that {E,g) itself, (and not just its tangent cones), is asymptotic 
to a fiat quotient of M.^, and hence has at most quadratic volume growth of geodesic balls or linear 
area growth of geodesic spheres. This is done in the following result, which is a strengthening of 
Lemma 3.2. 

Lemma 3.7. Under the assumptions (3.26) and (3.27) above, there is a compact set K C E such 
that {E \ K,g) is quasi-isometric to a flat product {M? \B) x or M+ x x S^. 

Proof: As in Case A, it is useful to work with the metric g = v?g\ again, this makes no 
significant difference, since (3.17) holds. All the metric quantities below are thus in the g metric. 
For notational simplicity however, we drop the tilde from the notation. 

Let t^{x) = lims^oo{dist{x, Se{s)) — s). As in the construction of Busemann functions, the 
limit here exists, c.f. [Wu] for a discussion of such functions. By construction, t^ is a Lipschitz 
distance function, i.e. too realizes everywhere the distance between its level sets. Observe that on 
M^,too is just the distance function to {0} e M.^ , on V x ,too is the distance function to {0} G V 
pulled back to y x 5^, for any cone V with vertex {0}, while on x x S^,too is the distance 
function on IR+ pulled back to the total space. 

By renormalization, (as with the potential u) , too induces a distance function too on each tangent 
cone ^oo by defining tooix) = lim{tooix)/too{xi)), where t{xi) oo and Xi are the base points 
converging to the base point Xqo G A^q. Thus too is the function above on x S*^ or R"'' x x S^. 

The map too '■ {E, g) — > M+ is distance non-increasing, and preserves distance along the integral 
curves of too; thus where smooth, it is a Riemannian submersion. We will show that too gives rise 
to a Lipschitz quasi-isometry by examining the asympotics of the second fundamental form of its 
level sets. 

Thus, let a{s) be any geodesic ray in E which is an integral curve of Vtoo, and let B = B^{s) 
denote the second fundamental form of the level surface to^^(s) at cf{s). The form B is well-defined 
and smooth along any such ray a. Recall that B satisfies the Riccati equation 

B' + B^ + Rt = 0, (3.28) 
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where T is the unit tangent vector along a. Consider the behavior of s-B{s) as s — > oo. This quantity 
is scale-invariant, and thus converges smoothly, (in subsequences), to the limiting expression s ■ 
Boo{s) on any tangent cone at infinity. Since the parameters s and too|<T are the same up to 
additive constants, s = too- Similarly, by the definition of too, -Boo is the second fundamental form 

of the levels too- Hence, either 

(i) : s- -Boo(s) = {d6/\d6\)'^, when the tangent cone is of the form V x S^, and 9 is the angle variable 
about {0} e V, or 

(ii) : s ■ Boo{s) = 0, when the tangent cone is of the form M"*" x x S^. 

Thus s ■ i?oo(s) is either of rank 1, with eigenvalue 1, or identically 0. Note that the expression 
in case (i) is independent of the cone V, i.e. the cone angle at {0}. 

As noted preceding Lemma 3.2 and in Remark 3.3(ii), the geometry of A{^s,2s) smoothly 
approximates that of a limit tangent cone Aoo, for s large, and varies continuously in s. Since the 
two alternatives (i) and (ii) above for the structure of B^o are rigid, it follows that all tangent cones 
^oo are of the same type, i.e. they are all of the form V X S^, or all of the form M"*" x x S^. 

The main task now is to show that the deviation of s ■ B{s) from its limit s ■ B^{s) has bounded 
integral. To do this, we use the Riccati equation, and estimate the decay of the curvature term Rt, 
using basically standard methods in comparison geometry, c.f. [P, Ch.6.2] for instance. 

Thus, from (1.4), the sectional curvature K of {M,g) satisfies 

i^xz = |Vi^P > 0, ifxy = -|Vi/p < 0, 

where Z = Vu/\Vu\, X,Y are vectors orthogonal to Z, and u = log u. Hence \Rt\ < |Vi^p. 
Substituting this in (3.28) gives 

\B' + B'^l <\Vi^f. (3.29) 

Let A be any eigenvalue of B, with unit eigenvector e; (note that B is symmetric). Observe then 
that s ■ A(s) converges either to 1 or to 0, as s — > oo. The estimate (3.29) when applied to (e, e) 
gives 

±(sA' + sA2) < s\Viy\'\ 
Integrate this by parts along any finite interval / to obtain 



|sA|a/ + J (sA^ - X)ds\ < J s\Viy\^ds. 

If sA ^ as s — > oo, choose / to be any interval on which sA = at dl and sA 7^ on /, so that s\ 
has a definite sign on I. If there are no such boundary points, choose I to be an infinite half-line. 
Similarly, if sX 1, choose / to be intervals such that sA = 1 at dl with sA — 1 7^ on /. Then 
summing up the estimate above over all such intervals gives 

/ \sX'^ - X\ds < / s\Vu\'^ds + Co, 

J (T J U 

for some constant Cg < 00. Now the estimate (3.18), together with (3.27) gives 

s\Vu\'^ds < Ci, 



for some constant Ci < 00. 

Suppose first s ■ X{s) — > 1 as s — >■ 00. We then obtain 



/ 



sA(A — -)\ds < C2 < 00, 
s 



and hence 



/ 



A - < C3 < 00. (3.30) 
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Similarly, if s ■ A(s) ^ as s — oo, one obtains 

j \X\ds < < oo. (3.31) 

Now the second fundamental form B gives the logarithmic derivative of the norm of Jacobi fields 
formed by the family of ioo-rays in E starting at some level t^{so)- Thus, if J is any Jacobi field 
formed from the too-conguence, and v = J/\J\, we have along any too-ray, 

B{v,v) = ^^{log\J\{s)). 

Suppose first that the end E has tangent cones at infinity of the form M"^ x x S^. Then (3.31) 
implies the uniform bound 

< \J\{s) < C4, (3.32) 

with C4 = This means that the geometry of the level surfaces of t^o is uniformly bounded as 
s 00, i.e. the diameter and area of the level surfaces is uniformly bounded away from and 00. 
It is then clear that there is a Lipschitz quasi-isometry of {E \ K, g) to M"'" x x induced by 
too- Since g and g are also quasi-isometric, by (3.17), the lemma follows in this case. 

If E has tangent cones at infinity of the form V X S^, then there is a basis of Jacobi fields whose 
elements satisfy either (3.32), or, from (3.30), 

Cl^s < \ J\{s) < C4S. (3.33) 

As before, this implies that too gives rise to a quasi-isometry of {E \ K,g) to (M^ \B) x S^. 

U 

Of course Lemma 3.7, in both cases, immediately implies that 

areaSE^s) < c - s, (3.34) 

for some constant c < 00. However, (3.34) violates the standing assumption (3.15), (or (3.2)). It 
follows that no end {E,g) can satisfy the assumptions of Case B. 

Together with Case A, this completes the proof of the second statement in Theorem 0.3. 

We now turn to the proof of the last statement in Theorem 0.3. We will assume that the end E 
is small, i.e. 



Ie 



-ds = 00, (3.35) 



Ie areaSE{s) 

and derive a contradiction from the assumptions supu < 00 and uie ^ 0- 

The proof is based on a result of Varopoulos [V] which states that ends of Riemannian manifolds 
satisfying (3.35) are parabolic, i.e. admit no non-constant positive superharmonic functions v which 
tend uniformly toward their infimum at infinity. (Actually, the result in [V] is a condition on the 
volume growth of geodesic balls, but this is equivalent to the bound (3.35) under the estimate 
(3.13)). We will prove that the potential u is such a non-constant function, giving the required 
contradiction. 

Thus, suppose first that 

supEU < 00, (3.36) 
Arguing as in (3.9), but now on £^ C M in place of N, we have 

|Vm|2 = / / \Vu\dA^dv, (3.37) 

E Jv JLyHE 
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and as in (3.10), 

/ |Vu| < f \Vu\ < C. 



iLvHE 

Thus, these estimates imply that 

/ iVup < oo, (3.38) 

J E 

SO that 

as s — > CO. (3.39) 

J Ae{s,x>) 

On the other hand, again referring to the proof of Lemma 3.6, since almost all (in terms of 
measure) points in n E may be joined by flow lines of Vu to a fixed bounded surface T in K ^ 
by the divergence theorem there is a subsurface T' CT such that 

[ |Vu| > [ |Vu| > c, (3.40) 

JL-^nE JT' 

for some constant c > 0. It follows from (3.37)-(3.40) that 

oscae{s,oo)U -^0 as s ^ oo. (3.41) 

By assumption (ii) in Theorem 0.3, we may thus assume w.l.o.g. that 

limt^ooU = 1, (3.42) 

in E. Thus, as noted in (3.18), |r| < e{t)/t^, where e{t) — >■ as t ^ oo, and so Lemma 3.2 holds 
on E. 

Now choose a smooth approximation S" to a large geodesic sphere Se{s) with Se{s) H dE = 0, 
so that 5* separates E into two components, one being the outside containing the end E. 
We claim that if, in addition to (3.36), 

m.E + 0, (3.43) 

then there is a set of flow lines 7 of Vtt or — Vu, starting on 5, of positive measure on 5, and 
pointing out of 5, which never intersects S again at later times. Hence such 7 diverge to infinity in 

since up to a set of measure 0, 7 does not terminate in a critical point of u. (Compare with the 
earlier discussion regarding divergent flow lines and Weyl dipole solutions concerning Lemma 3.4). 

To see this, suppose instead that all flow lines say of Vu which initially point out of S eventually 
intersect S again, with the exception of those terminating in critical points. Consider the measure 

=< Vu, V > dA, (3.44) 

on S, where dA is the Lebesgue measure and u is the unit outward normal on S; dji is absolutely 
continuous w.r.t. dA. Since u is harmonic, the divergence theorem implies that the gradient flow of 
u preserves the measure d// in the following sense. Let D be a domain in S and let O be the domain 
in E formed by a collection of flow lines outside 5, whose endpoints form a smooth surface D' . If 
v' denotes the unit outward normal to Q, at -D', then the flow from D to D' carries the measure dji 
to the measure dji' =< Vu, v' > dA', where dA' is Lebesgue measure on D'. In particular, the flow 
preserves the masses of the measures. 

Thus, under the assumptions above, the gradient flow (with varying flow-times), induces a home- 
omorphism of S'\ Z into itself, where Z is a set of Lebesgue measure 0, corresponding to flow lines 
terminating in critical points. However, this homeomorphism inverts the direction of Vu w.r.t. the 
fixed normal v, and hence maps domains on which the measure d/x is positive onto domains 
-D- on which d/j, is negative, in such a way that that 

m^(L>+) = |m^(Z?_)|. (3.45) 
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This of course implies that the total mass of d/j, on S is 0. However using (3.44) and the remarks 
concerning (3.25), the mass niE of E equals the total mass of d/i. This contradiction proves the 
claim. 

We may now complete the proof as follows. Assuming E is an end satisfying (3.36) and (3.43), 
there exists an open set O of flow lines 7 = 7(t) of either V?/ or — Vu which start on a set of 
positive measure on S and diverge to infinity in E. Consider the former case, which corresponds 
to > 0. A generic flow line of Vu in E tends to the maximal value of u in E and hence a 
generic flow line in O also tends this maximal value. By (3.42), it then follows that supgy/^u = 1, 
for some compact set K (Z E. The function v = —u is thus a bounded (non-constant) harmonic 
function on E, which tends uniformly to its infimum at infinity. Hence E cannot be parabolic. 
This contradiction shows that (3.35) cannot hold for E^ and thus, by the proofs in Cases A and B 
above, the end E is asymptotically flat. The proof in case ruE < is the same. 

This completes the proof of Theorem 0.3. ■ 

Remark 3.8.(1). There are (non-flat) static vacuum solutions with a small end, namely the 
Kasner metrics (2.33), with a > 0. These solutions have volB(s) ~ s'^~^, volS{s) ~ s^~^, where 
(5 = (a + — 1)"^ G (0, 1), and hence the end is small. Note that u ~ s'^ is unbounded. There 
are other Weyl solutions which are complete away from E with dM pseudo-compact and with one 
small end, (take for instance the potential given by the Green's function on x S"-*^, see (ii) below), 
but all known examples with small ends are either asymptotic to the Kasner metric at infinity or 
have faster than quadratic curvature decay, i.e. satisfy (3.18). 

It is an open problem to understand in more detail the structure of small ends of static vacuum 
solutions. It follows from the results above in §3 that all tangent cones at infinity of E are collapsing, 
and hence they are all Weyl solutions. But the metric uniqueness of tangent cones at infinity is 
unknown, as is the question of whether small ends have finite topological type. 

(ii). We construct an example which illustrates the sharpness of the last statement of Theorem 
0.3. Let Gi and G2 be the Green's functions for the Laplacian on the flat product x S*^, with 
poles at {0,pi), for pi,p2 distinct points in S^. Here we consider as the z-axis in quotiented 
out by an isometric Z-action. As in (2.34), Gi{x) = G{x,pi), viewed as a function on the universal 
cover R^, may be written as 

n ^ 

Gi = limn-,00 {SA6) 

j=-n -I 

where rj is the Euclidean distance to the collection of lifts of pi in W"' and c„ is a suitable normalizing 
constant with c„ 00 as n — >■ 00, chosen so that G{x,pi) is finite. Thus Gi is an axisymmetric and 
z-periodic harmonic function on (R^)+ x S^, where the now means rotations about the z-axis 
in R^. Hence the potential 

u = Gi-G2 

generates a Weyl solution as in (2.12), with u = and which has an isometric Z-action along 
the 2;-axis. Let M denote the Z quotient of this solution. Then the metric boundary of M is 
pseudo-compact. Since Gi ~ logr as r — > 00, 

n — > 1, 

at infinity in M. The end E = R"*" x is small and has mass in the sense of (0.7). These 
solutions of course resemble the dipole-type solutions discussed in §2(IIC), but with a collapsing 
end. 

Remark 3.9. (i). Although we will not detail it here, an examination of the proof shows that 
Theorem 0.3 holds for non- vacuum static solutions of the Einstein field equations (1.2), provided 
suitable decay conditions are imposed on the stress-energy tensor T. 
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This is the case for example, if T is any tensor with compact support, or more generally if T 
satisfies an estimate of the form |T|(x) < c-t^^(x), for some c > and all x with t{x) > Sq, together 
with ^Au G L^{M \ U). This latter condition is needed to obtain the bound (3.10). By (1.2), note 
that ^An = ^trT. The starting estimate (3.3) may be obtained by applying (a suitable version of) 
[Anl,Thm.3.3]. 

(ii). Also, Theorem 0.3 can be given a finite or quantitative formulation, i.e. one can relax the 
assumption of completeness, in basically the same way as the local estimates (3.3) follow from the 
non-existence of global static vacuum solutions with dM = (Ji, c.f. [Anl.App.]. 

Thus, if (3.15) holds and {M,g) is 'sufficiently large', depending on K, then sufficiently far out 
in {M,g), the metric is close to a flat metric. We leave a precise formulation and proof, (based on 
Theorem 0.3), to the reader. 

Remark 3.10. As noted in §0, it would be of interest to prove that {M,g,u) has a unique end. 
Under the hypotheses of Theorem 0.3, wc conjecture this is the case at least when M is complete 
away from S. If M is in addition smooth up to S, this has been proved by Galloway [G]. Following 
essentially the same arguments as in [G] , it is not difficult to show that if M is complete away from 
S and if the Riemannian 4-manifold N admits a compact smoothing of a neighborhood of S having 
non-negative Ricci curvature, then M has a unique end. 

Remark 3.11. We point out that Theorem 0.3 and Theorem 0.1 are false in higher dimensions, 
due to the existence of Einstein metrics on compact manifolds, which arc not of constant curvature, 
in dimensions > 3. (The equations (0.1) on any n-dimensional manifold Af " generate Ricci-flat 
manifolds on N'^'^^). 

Thus, let (E, g) be any compact (n — 2) dimensional Einstein manifold, with Ricg = {n — 3) ■ g 
and define the warped product metric ^ on x S by 

[n — ly 

where / is the unique function on [0, oo) such that /(O) = 1, /' > and (/')^ = 1 — /^~"- A simple 
computation shows that (M^ x S,^) is complete and Ricci-flat, c.f. [Bes, p. 271] and the space-like 
hypersurface M"*" x S, with metric dt^ + /^(t) • 5, is a solution to the static vacuum equations, 
with potential u = /', (up to a multiplicative constant). Thus the horizon is S and the solution is 
smooth up to E and complete away from S. 

This metric is asymptotically conical, i.e. asymptotic to the complete (Euclidean) cone on (S,^), 
but is asymptotically flat only in the case that (S, 5) = ^""^(l), corresponding to the n-dimensional 
Schwarzschild metric. 
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